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The  theory  of  the  computable  a-posteriori  error  estimate  for  a  finite  element 
method  is  developed.  Among  other  things,  it  is  shown  that  the  error  estimate  is 
very  reliable  and  the  ratio  (called  effectivity  index)  between  the  estimator  and 
the  true  error  approaches  one.  Numerical  examples  computed  by  program  FEARS 
(Finite  Element  Adaptive  Research  Solver)  of  the  University  of  Maryland,  illustrate 
the  effectivity  and  reliability  of  the  estimators. 


1.  INTRODUCTION 


Recent  Iv  an  increasing  interest  in  the-  finite  element  computations  is 
being  focused  on  the  reliability  of  the  results  and  the  quality  of  the  used 
meshes  and  elements. 

During  recent  years  at  the  University  of  Maryland,  the  studies  were  under¬ 
taken  which  focused  toward  the  development  of  a  finite  element  system  having 
the  following  features. 

a)  The  solver  supplies  the  user  with  a  reliable  and  accurate  information 
about  achieved  accuracy  in  the  desired  norm. 

b)  The  solver  constructs  adaptively  meshes  which  are  leading  to  the  highest 
possible  accuracy  (through  an  adaptive  refinement). 

c)  The  solver  uses  the  most  simple  input. 

d)  The  solver  combines  the  advances  in  the  mathematics  and  computer  science 
including  parallel  computations. 

The  solver  FEARS  (Finite  Element  Adaptive  Research  Solver),  its  mathematical 

version  FM  developed  for  Univac  Series  1100  implements  some  of  the  points  mentioned 
* 

above  .  The  detailed  description  of  FEARS  and  the  experience  with  it  will  be 
published  elsewhere.  For  some  information  about  FEARS  and  its  applications, 
we  refer  to  |  ;>  ]  ,  f  3  ]  ,  |  4  ]  ,  [  5  1 . 
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For  analysis  of  the  pnrnllelity  we  reter  to  [  1  ]. 
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One  of  the  main  aspects  of  FEARS  is  the  theory  of  the  a-posteriori  estimates. 

Some  aspects  of  the  a-posteriori  error  estimates  and  optimality  of  the  meshes 

were  investigated  in  [6  1  ,  [  7  1  ,  1  8  ]  . 

The  error  norm  | |a||  is  approximated  by  the  computable  estimator  e 

computed  through  error  indicators  n^(A)  associated  to  every  element  A  and 

computable  locally  by  knowledge  of  the  finite  element  solution  at  the  particular 

£ 

element  A  and  its  direct  neighbors.  The  effectivitv  index  0  =  | | — |  j-  expresses 
the  quality  of  the  estimator  and  0  should  be  close  to  one  when  the  error  is 
sufficiently  small  (e.g.  5/1).  It  is  desirable  that  the  estimator  z  has  the 
following  two  properties: 

(1.1)  0  *  CL  £  6  1  Cy  <  » 

with  C  and  C.,  independent  of  the  solution  and  the  meshes  under  very  general 
condit ions . 

(1.2)  0  ->  1  as  |  |e|  |  -*•  0 

provided  that  some  additional  assumptions  about  smoothness  are  made.  The  present 
paper  develops  the  theory  of  the  estimator  which  satisfies  (1.1)  and  (1.2),  and 
is  implemented  in  FEARS.  The  energy  error  norm  is  assumed  and  model  elasticity 
problem  is  considered. 

Section  2  consists  of  some  preliminary  notions. 

Sections  3  and  4  elaborate  on  the  type  of  meshes  which  are  adaptively 
const  meted . 

Section  5  deals  with  the  approximation  properties  of  the  elements  on  the 
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admissible  meshes. 

Section  6  formulates  the  model  problem  (elasticity  problem). 

Section  7  develops  the  estimator  and  proves  (1.1). 

Section  8  proves  that  the  estimator  is  asymptotically  correct,  i.e. 

0  ->  1  . 

Section  9  deals  with  two  computational  examples  and  discusses  the  effectivity 
of  the  approach. 

The  adaptive  construction  of  the  meshes  is  based  on  the  equilibration  of 
the  error  indicators.  This  principle  was  theoretically  analyzed  in  [  7  ]  for 
one  dimensional  problems  and  its  theoretical  investigation  in  the  context  of 
FKARS  will  appear  elsewhere. 


I 

i 
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2.  BASIC  NOTATION 


2 

Throughout  this  article  we  denote  by  R  the  two  dimensional  real  Euclidian 

2  2  2^/2 

space  with  x?  (x-pX^GR  ,  ||x||  =  max(  j  x1 1  ,  |  x2  |  ),|  |  x  |  |  E  =  [x^+x^]  .  Let 

2 

QCR  be  a  bounded  set  and  3Q  its  boundary.  We  define 


diam  Q  =  sup  |  j  x-y |  | 
x,v6Q 


dist(x,Q)  =  inf | | x-y | | 
v€Q 


and  for 


Q.€  R 


1,2 


dist(Q  ,Q„)  =  inf  |  |  x-y  |  | 
x€Q2 


E 


An  index  E  will  denote  that  the  norm  ||'||  is  used  instead  of  ||‘|| 

,  dist  (x,Q)  =  in f | | x— v |  | 

y€Q  E 

For  a  p  >  0  ,  QC  is  the  p-neighborhood  of  Q 


0 


P 


{ x€TR2 


dist (x,0) 


<  P 


}  . 


The  closure  of  Q  in  IR2  is  denoted  by  Q  ,  int  0  means  as  usual  the  interior 
of  Q  . 

2 

By  7.  we  denote  the  set  of  all  two  dimensional  integers  k  i  (k^,k2>  , 
kj ,  i  =  1 , 2  integral . 

2 

Suppose  6  >  0  is  a  positive  real  number,  then  we  will  write  for  any  k€Z 
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Qk  =  {x€IR2 Ik. 0  <  x.  <  (k.+l)6  ,  i  =  1,2} 

e  1  i  -  i  -  i 

2 

Assume  that  Z  C  Z^  is  a  finite  set.  Then  we  denote 
o 


int  [ 


U  Qgl 

k€Z  9 
o 


We  shall  assume  that  is  such  that  0  is  a  Lipschitz  domain.  For  brevity, 

o’ 

whenever  it  cannot  lead  to  misunderstanding  we  shall  write  instead  of  iL  q  . 

2  ° 

When  we  talk  of  a  square  S  in  IR  ,  we  shall  always  suppose  that  it  is  closed 

and  that  its  sides  are  parallel  to  the  coordinate  axes,  i.e.  S  is  of  the  form 

[a,a+dl  x  [b,b+d]  for  a,b,d£IR,  d  >  0  . 

As  usual,  let  L?(Q)  =  H°(f2)  be  the  space  of  all  square  integrable  functions 

on  SI  with  the  inner  product 


(u,v) 


L2(fi) 


uvdx,  dx  =  dx^dx2 


and  the  corresponding  norm  | ! 
the  usual  Sobolev  space  with  the  norm 


l2(w)  . 


By  H  (W)  ,  k  >  0  integral  we  denote 


H  (Q)  0 <  |  <x  |  <k 


.1.  llDau|| 


L2(n)  * 


where  x  =  (a^,«2),  a.  >  0  ,  |a|  =  +  ot 2 


and 


pa  = 


ai  11 2 
9x^  3xz 
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Obviously  we  have  H°(ft)  =  L2(ft)  .  We  will  also  use  the  notation  j 

i 

s 

t 

I 

i 

f 

in  the  usual  (distributive)  way  ' 

be  the  completion  of  the  set  of 

We  will  also  deal  with  functions  defined  on  one  dimensional  manifolds, 
more  precisely  on  the  boundary  3ft  or  a  part  T  of  it.  The  notation 

L?(F)  =  H°(r)  has  then  the  obvious  meaning. 

m  k(i) 

Let  F  =  .U  r  ,  where  each  F.  is  a  closed  side  of  some  0A  Cft  , 

1=1  1  l  <6 

with  r.C3^  ,  (i  =  l,...,m);  then  we  shall  write 
ll1,r(Sl)  =  MiCH1  (ft)  u  =  0  on  n  . 

Obviously  (ft)  -=  H1,r(ft)  when  F  =  3ft  and  H1(fi)  =  H1,F(ft)  when  T  =  0  . 

Finally  by  C°(ft)  we  denote  the  space  of  all  continuous  functions  on  ft 
and  let 

I  |  u  |  |  =  sup | u (x) I 

C°(ft)  x€ft 

1  0 1 

We  will  deal  later  with  extensions  of  functions  in  H  (ft),  H  (ft)  and 
H^’1 ( / )  from  into  a  neighborhood  of  ft  . 

Theorem  2.1.  There  exists  an  operator  T  mapping  (ft)  into  H^(ftP), 

(where  ft'  is  a  e-neighborhood  of  ft.  )  such  that 


l"l\  -l  IlnMI,2  (u) 

H  (ft)  |af=k  V  '' 


We  define  the  support  of  a  function  u€L2(ft) 

°k  k 

and  denote  it  by  supp  u  .  Let  H  (ft)cH  (ft) 


all  functions  having  compact  support  in  ft 


i)  for  a  2<6<°°,  0  <  <  p  and  any  square  SCfi  and  any 

0  <  a  <  (\  we  have 
o 


|Tu||  i  a  1  C | I u l i 

H1(SJl)  H1(S6ariSi) 


with  C  independent  of  a,  S,  u 

ii)  If  x€!2P  ,  x^ft  ,  T  ^  ,<ft 


and  distE(x,P)  £  Jq  dlstgCx, 9ft-f)  then 


Tu  =  0  on  (x)-ft,  y  =  dist_(x,r) 

Z  L 


with  S  (x)  being  the  square  with  the  center  in  x  and  lengthside  y 


then  Tu  =  0  on  ftP  -  ft 


iii )  If  I 

= 

9ft  i.e 

P_roof.  I f  is  enough  to  prove  the  theorem  in  the  neighborhood  of  the  boundary 

(i.e.,  endpoints)  of  F  .  In  the  neighborhood  of  all  other  points  x€9ft  we 

use  the  classical  extension  theorem  when  x^F  and  we  extend  by  zero  when 

x  €  1  and  apply  the  standard  argument  with  partition  of  unity. 

The  endpoint  of  T  can  be  located  in  a  vertex  of  9ft  with  the  internal 
3  1 

angle  ^  "  or  f  or  it  can  be  on  straight  part  of  the  boundary. 

We  will  deal  only  with  the  case  of  the  vertex  in  the  coordinates  origin 

3 

and  the  internal  angle  on  n  .  The  other  cases  are  analogous. 

Let  (see  Fig.  2.1) 


ftP  =  (x€ftf  |  x^ft,  dist  (x,9ft-D  _<  4  dist„(x,  T)  ) 


sclirmi 


notation  in  the  proof  of  the  Theorem  2.1 


<a 
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Obviously  ^***  consists  of  a  sector  with  the  lines  and  as 

its  boundary.  Let  V  be  a  symmetric  (with  respect  to  the  origin)  sector 

%  'v,  i); 

with  the  boundary  S1  and  and  let  be  the  line  symmetric  to 

.  % 

(with  respect  to  ?<?.  -  V)  .  Finally  let  VCS1  (respectively  VC2)  be  the  sector 

bounded  by  312  -  T  and  (respectively  and  W  be  the  sector  bounded 

>yl 

by  S  and  F  . 

1  r 

Assume  now  that  u  €H  ’  (12)  .  By  an  afine  transformation,  we  construct 
w£H  (V)  such  that  w  =  u  on  3£2  -  T  ,  w(x)  =  u(x)  for  x€S^,  x€S^,  and 
j|x||p  =  ]  I x I  I E  .  Obviously 


3) 


l1  (.V) 


The  extention  Tu  on  be  now  the  reflection  of  w  (with  respect  to 

■K:  -T)  .  It  is  readilv  seen  that 


i  I  Tu , 


V(* 


1  1  ^ 
HJ(V) 


and  f or 

x€St,  x^,  I  I  Jt !  |£  =  I  |x|  |E 

we  have 


(Tu) (x)  =  u(x) 


We  extend  u  on  by  zero. 

Bv  an  afine  transformation  of  W  on  V  we  can  easily  continue  a  function 
v  such  that  v  =  u  on  and  v  =  0  on  and 
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H1  (V) 


! u  i 


H  (W) 


Now  let  Tu  on  1  be  the  symmetric  image  of  v  .  It  is  easy  to  see 

that  our  construction  has  all  properties  of  the  extension  formulated  in  the 
theorem  when  is  chosen  sutticientlv  large. 


► 

:  j 

'  I 
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3.  THE  MESH  AND  ITS  BASIC  PROPERTIES 

Wo  will  introduce  now  a  class  of  partitions  of  SI 

'o  ’ 

We  define  a  mesh  V(i!)  =  {A  }  as  a  finite  collection  of  closed  squares 
AlCfi  of  various  sizes  with  sides  parallel  to  the  coordinate  axes,  and  which 
are  generated  by  the  following  recursive  rules. 

k 

i)  The  squares  (.Q^l  >  k €  Z  ,  create  a  mesh. 

ii)  If  {A1}  ,  i  =],..., m  is  a  mesh,  then  a  new  mesh  is  obtained  if  any 
A1  is  subdivided  into  four  congruent  squares  of  half  the  side  length  of  A*  . 

Any  &X€D  will  be  called  an  element  and  its  sides  the  edges.  The 
vertices  of  the  elements  will  be  called  the  nodes.  A  node  P  will  be  called  a 
regular  node  if  either  P6  3.Q  or  P  is  a  vertex  of  four  different  elements. 
Otherwise  P  is  an  irregular  mode.  By  P( V)  we  denote  the  set  of  all  nodes 
of  V  and  by  R(V) C  P(V)  the  set  of  all  regular  nodes.  Finally  let  h(P)  = 

ma''  diam  . 

A€P 

Figure  3.1  shows  an  example  of  a  mesh.  The  irregular  nodes  are  marked  by 

a  cross. 


Figure  3.1.  An  example  of  a  mesh. 
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Let  D  ’.A1)  ,  i  =  Then  obviously  Jj  A1  =  U  .  We  shall 

i  =  l 

denote  by  AI(D)  ,  the  subspace  of  all  continuous  functions  on  fl  which  are 
Individually  bilinear  on  each  A*  ,  i  =  1 , . . . ,m  .  It  is  clear  that 
M(P)CH1r;.)  • 

We  will  always  assume  that  at  least  four  different  h€V  lie  in  every 

Q-  • 

Now  we  will  analyze  further  the  basic  properties  of  the  meshes  introduced 
above . 

LEMMA  3.1  .  Let  A  '  ,A"€  V 

1)  Assume  that  £  1  f)  A"  4-  0  and  diam  A'  <-  diam  A"  .  Then  one  and  only 

1  'nc'  of  the  following  state  men  t  s  holds  . 

i) 

i 1 )  A 1  fl  "  is  just  one  point  being  a  common  vertex  of  A’  and  A"  . 

i ii)  A  1  fl  A"  i s_ _a n_  edge  of  A 1  and  is  contained  in  an  edge,  T  say. 

of  A"  .  If  x' ,  x"  ere  any  two  endpoints  of  (A'PlA")  and  F  respectively, 

then  |x'-x”|  is  an  integral  multiple  of  diam  A'  . 

2)  If  P  €P(P)  ,  P  ^  R  (V)  (i.o.,  P  is  an  irregular  mode)  and  P  is  a 

vertex  of  A'  €  V,  then  there  exists  A"€  D  such  that 

i  )  P €  A" 

Li)  diam  A"  diam  A' 
iii)  A”  0A'  is  an  edge  of  £.'  . 

3)  diamA"/diamA '  =  2s  with  s  an  integer. 

The  lemma  can  be  ..asily  proven  by  induction. 

LEMMA  3.2.  Suppot  e  that  AJz  D^. _ Then  at  least  one  vertex  of  A  is  a 


regular  node  and  If  ari  edge  of  A  is  contained  in  30  for  some  k€Z  , 

_ _ _ _ _ _ _ _ _ _  'Q _ o 

k 

then  at  least  one  vertex  of  A  which  is  not  on  is  a  regular  node. 

Proof .  Observe  that  any  regular  node  always  remains  regular  when  our 
recurrent  construction  is  implemented.  The  lemma  follows  now  easily  since  at 
each  step  the  midpoint  of  the  subdivided  element  becomes  a  regular  node  and  it  is 
a  vertex  of  all  four  new  elements  created  at  that  step. 

As  seen  in  Figure  3.1  there  could  be  an  element  A€P  ,  such  that  only  one 
of  its  four  vertices  is  a  regular  node. 

LEMMA  3.3.  Let  V  be  a  mesh  and  P€R(P)  .  Then  there  exists 
Vp€M(P)  such  that  v^  (P)  =  1  and  v^CQ)  =  0  for  any  Q€R(fl)  ,  Q  4  P  . 


Proof.  We  first  note  that  it  suffices  to  define  v  at  the  vertex  of  each 

-  p 

&€V  .  Let  V  =  {A1}  ,  i  =  l,...,m  .  Assume  that  we  have  denumerated  the 
elements  so  that  diam  A1  >_  diam  Ai+1  i  =  l,...,m-l.  We  prove  now  our  lemma 
by  induction  with  respect  to  i  . 

First  let  us  observe  that  all  four  vertices  of  A1  are  regular  nodes. 

If  one  were  irregular  then  by  lemma  2.1,  there  exists  A€P  such  that 

%  l  1 

diam  A  >  diam  A  .  This  is  a  contradiction  because  A  is  the  largest  element. 

Suppose  now  that  j  =  1  ,  or  j  >  1  and  v  has  already  been  constructed 
i= j -1  .  .P 

1  1  T 

on  U  A  Consider  now  the  vertices  of  A'  .  If  Q  is  a  vertex  of  AJ 

i=l 

and 

i)  Q  is  a  regular  node,  we  define  v^CQ)  =  for  Q  =  P  and  v^CQ)  =  0 
for  Q  /  P  . 

ii)  Q  is  an  irregular  node,  then  by  Lemma  2.1,  Q€A'  where 

diam  A'  diam  A*  .  By  induction  assumption  v^  was  already  defined 

on  A'  and  so  v  (Q)  is  defined  too.  We  therefore  construct  v 

i  =  l  ^  p 

on  A^  with  desired  property  on  U  A  Let  us  remark  that 

i=l 
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i=M  i 

if  o£  U  a  then  the  "a  I  ue  v  (0)  =  1  respectively  0  is  the  same  as  in 
i  =  i  p  ' 

the  previous  phase. 

LEMMA  3.4.  Let  D  he  a  mesh  and  u£M(P)  be  such  that  u(Q)  =  0  for  any 
Q€  R(V)  then_  ju  =  JJ_  . 

Proof.  Let  { A 1 i  be  numerated  as  in  the  proof  of  Lemma  3.3.  The  lemma 

will  be  proven  by  induction  with  respect  to  i  .  Because  all  vertices  of  A^ 

are  regular  nodes  we  have  u  =  0  on  A ^ .  Let  now  i  =  l,...,j-l  <  m  and  consider 

u  on  A-*  .  if  the  vertex  P  of  A~*  is  a  regular  node  then  by  assumption 

u(P)  =  0  .  If  P  is  irregular  node  then  by  lemma  2.1,  we  have  also  P  €.  A1 

with  diam  A*  diam  .‘r  ,  i.e..  A1  =  A^  ,  for  some  k  '■  j  .  By  induction 

1 = i  - 1  i 

u  =  0  on  .U,  A  and  therefore  u(P)  =  0  .  So  u  =  0  in  all  vertices 
i  =  l 

of  and  so  u  =  0  on  A*  and  lemma  is  proven. 

Lemmas  3.3  and  3.4  show  that  the  function  n€M(P)  is  uniquely  defined 
hv  its  values  at  t  lie  regular  nodal  points. 

Lemma  I. l  and  arguments  analogous  to  those  used  in  the  proof  of  Lemma 
3.3  yield 

LEMMA  1.3.  Let  l'  he  a  mesh .  Then 

i)  The  set  of  functions  1 v  | P € R (V) (  creates  a  basis  for  M (V) 

i  i  )  v  ■  >') 

P 

i  i  i )  )  v  =  1 

P  €  R(P)  r 

Definition  3.J:  The  set  =  supp  ,  J’€R(P )  will  be  called  a  star 

associated  to  the  node  P  or  briefly  a  P-star. 

_  f  1 

P 


lemma  3.3  (iii)  yields  readily  that  U 

P  €R(V) 
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LEMMA  3.6. 

Let 

A  €V  , 

P  €  R(0)  and 

A  fl  int  o! 

_  _  P_ 

i  0  . 

Then  ACu  . 

P 

Proof.  Assume 

on  the  contrary  that  A 

<£(o  .  Then  there 

~  P 

exists  an  open 

set  SCA  ,  S  fl 

P 

0  .  Because  v  £M(D) 
P 

,  v  is 

P 

bilinear  on  A  .  But 

v  =0  on  S  ,  and 

P 

hence  v 

P 

=0  on  A  . 

This  is  a 

contradiction  because 

we  assumed  that 

A  fl  int  w  ^ 
P 

0  . 

LEMMA  3.7. 

The 

se  t  m 

P 

is  connected 

in  the  sense  that 

for  any  two  elements 

A',  A"C  ui  there 

P 

exists  a  sequence  of  elements  A'  = 

&0’A1” 

.  .  A  =  A"  such 
n 

that 

i)  A . C  ui 

1  P 

ii)  A  .  fl  A.  , 
l  l+l 

*  0 

and  is 

the  edge  either  A^  or 

Ai+1 

Proof.  Let 

OJ 

S  i 

=  u  A-1 

,  A 1  £  V  .  Assume  that 

we  have 

enumerated  the 

elements  A*  of 

P 

'0 

P 

so  that 

diam  A-'  >  diam  A^  +  1  , 

.1  =  1  >  • 

. . ,s-j  .  It  is 

'  1 

obvious  that  we  can  restrict  ourself  to  the  case  when  A  =  A 

First  we  prove  that  one  of  the  vertices  of  is  the  node  P  .  Suppose 

that  the  node  P  is  not  a  vertex  of  .  Let  Q  be  any  vertex  of  A^  .  If 

Q  is  regular,  then  v  (Q)  =  0  .  If  Q  is  irregular,  then  there  exists 

f\*€V  with  Q  £  A*  and  diam  A*  >  diam  A^  .  Thus  A*  flint  w  =  <J)  and  so  v  =  0 

P  P 

on  A*  by  lemma  3.6  and  it  follows  that  v^(Q)  =  0  ,  so  Vp  =  0  on  A^ 

which  is  a  contradiction.  Now  we  prove  the  lemma  by  introduction  with  respect 

k  1 

to  j  .  Assume  therefore  that  we  are  able  to  connect  A  with  A  and 

k+ 1  k+ 1 

consider  the  element  A  .  If  a  vertex  of  A  is  the  node  P  ,  then 

k+1  1 

A  can  be  connected  with  A  because  fi  is  a  Lipschitz  domain.  If  all 

k+1  k+lj. 

vertices  of  A  were  regular  nodes  different  from  P  then  A  <tw  by 

’  P 

lemma  3.6. 


So  we  need  onlv  consider  the  case  where  a  vertex  R  of  A 


k+1 


is  irregular 


and  v  (R)  ^  0  .  Bv  lemma  3.1  there  exists  A* €  V  such  that  diam  A*  > 

P 

k+1  _  k+1  W+1 

diam  A  ,  A*Oa  is  an  edge  of  A  and  R  €  A*  .  Thus  A*Cu  and 

P 

so  A*  =  a'  ,  j  '  k  .  So  A^+^  can  he  connected  with  A^  and  therefore 

with  A^  and  lemma  is  proven. 

Lemma  3.7  shows  that  int  w  is  a  domain. 

P 

So  far  we  have  not  made  any  restrictions  concerning  the  mesh  V  .  In 
the  next  section  we  will  analyze  the  family  of  K-meshes,  which  play  an 
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A .  THE  K-MESH 


Del  in  it  ion  A . 1 .  Lot  K  •  0  ,  real.  A  mesh  V  will  be  called  a  K-mesh 


if  for  any  I1  €  R(V) 


( A . 1 )  d i am  ,  <  K  inf  d i am  A ’ 

P  A  '€V 

A '  C  i. 


The  definition  has  a  clear  sense  because  of  Lemma  3.6. 

We  conjecture  that  definition  A.l  is  equivalent  to  3K*  >  0  such  that 

c  1 1  a  diam  A  , 

for  all  A  €  V  ,  sup  — - — .-  <  K* 

.  n  .  i  .  ,  c  »t  diam  A  — 

A  lift  is  an  edge  of  A 

Everywhere  in  what  will  follow  we  shall  assume  that  we  deal  only  with 
K-meshes.  We  mostly  will  not  mention  it  explicitly. 

LEMMA  A.l.  Suppose  V  is  a  K-mesh.  Then  there  exists  numbers  Mr  N 
depending  only  on  K  such  that 


i) _ I  f_  ?€R(P)  then  the  P-star  consists  at  most  of  N  different 

elements  of  V  . 

ii)  If  A'  €  V  then  A'Crn^  for  at  most  M  different  F€R(V)  . 

iiij _ I  f_  then  there  are  at  most  AK  +  A  elements  P  such 

that  A ' 0  A"  i  0  . 


Proof . 

i)  The  star  u<  can  be  contained  in  a  square  S  with  its  side  diam  w 
P  P  P 

So  for  the  number  N  of  elements  contained  in  oi  we  have  the  simple  estimate 

P  P 


area  w  fdiam  m  1 

_ _p  _ _ p 

inf  (area  A')  -  [inf  (diam  A ' ) ]^~  -- 

A  '€  '0  A 
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The  node  F  must  he  a  vertex  of  some  A" C  o>  satisfying  (4.2).  The  number 
of  such  elements  is  obviously  bounded  by 


area  Q 

_ EL-  < 

.diam  A '  — 

(  K  ' 


4K2(K-l/2) 2 


Because  not  more  than  4  regular  nodes  could  be  on  any  element,  we  see  that 


M  <  16K2(K-l/2)2 


iii)  Any  vertex  of  A'  can  be  a  vertex  of  at  most  three  other  elements 

and  therefore  by  lemma  3.1  it  is  sufficient  to  bound  the  number  of  elements 

A"  with  diam  A"  diam  A'  such  that  A'^flA'  is  an  edge  of  A"  contained 

in  some  edge  of  A '  ,  say  P  .  Fix  this  T  ,  and  suppose  that  there  are  q 

such  \"  .  q  is  finite  since  by  the  definition  of  V  ,  there  is  only  a 

% 

finite'  number  of  elements.  Then  for  at  least  one  such  A"  ,  A  say  ,  diam  _ 


—  diam  A' 

q 


Since  J  v  =  1  by  lemma  3.5,  there  exists  P€R(V)  such 
P€R(D)  p 


%  _  -v 

that  v  is  not  identically  zero  on  ALIA'  ,  and  therefore  both  A  and 
P 


A'C 


But 


v,  % 

q  diam  A  <  diam  A'  <  diam  to  <  K  diam  A 

p  _ 


Hence 


q  ;  K 

and  the  lemma  easily  fellows. 


LEMMA  4.2.  There  exists  a  number  L  (depending  only  on  K)  such  that 


OP 
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R(P)  can  he  partitioned  into  j  <  1,  sets  x  j  >  j  =  such  that  if 


P,o€  v 


J 


^  Q  then  int  ui^Dint  i,|q  =  0 


Proof.  Suppose  P  4  Q  ,  P,()€R(D)  and  int  u)p  flint  Wq  4  0  •  By  lemma 

3.6  int  t.)  flint  au  must  contain  the  entire  interior  of  at  least  one  element. 

P  Q 

Bv  Lemma  4.1  tli  contains  at  most  N  elements  and  each  of  these  elements 
P 

can  be  contained  in  o  ,  for  at  most  M-l  nodes  P' €  R(V)  ,  P  4  P'  •  There- 

P 

fore  there  can  be  at  most  N(H-l)  regular  nodes  Q  such  that  int  Wpflint 

We  shall  construct  now  the  sets  Xj  by  the  following  recursive  procedure. 

Let  P  ,  P?,..,Pr  be  some  enumeration  of  the  regular  nodes.  Set  =  {P^K 

l- l  l-l 

Suppose  that  we  have  already  defined  sets  , .  . .  ,  Xg  for  some  i  1-  » 


l-l 


l  •  1  1.  [f  for  P^,€R(P)  and  some  1  £  k  <_  S^_^  ,  int  u>  flint  0)^-0 


for 


all  Q  € 


l-l 


,  then  choose  k  to  be  minimal,  and  set 


/  =  J  for  t  €(],...  ,S^  -  (k) 


■i  ■  4'1|jiv 


and  define  =  S^_^  .  Otherwise  set 


l  £-1 

,t  =  for  1  <  t  <  S 


£-1 


>S£-1  +  1  =  {P£} 


w  i  tii  =  S^_  j  +  1  . 

Now  from  the  first  part  of  our  proof  we  see  that  N(M-l)  +  1  and  so 

L  <  N ( M- 1 )  +  1  and  the  lemma  is  proven. 
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Remark.  As  we  mentioned  in  the  introduction,  this  paper  develops  the 
basic  ideas  of  [1]  .  The  lemma  4.1  relates  to  the  intersection  index  and  the 
lemma  4.2  to  the  overlap  index  as  introduced  there. 

Given  an  element  A€P  we  will  always  enumerate  its  vertices  as  shown 
in  Figure  4.2. 


Figure  4.2.  The  Numeration  of  the  Vertices 


Let  P€R(D)  and  w  its  star.  The  node  P  is  a  vertex  of  at  most  4 

P 

elements.  We  will  denote  by  A  cm  the  element  for  which  the  vertex  number 

r  r 

given  by  Fig.  4.2  has  the  minimal  value.  This  rule  associates  to  every  P  a 
unique  element. 

Let  now  u  be  a  P-star.  Then  by  J  we  denote  the  invertible  affine 
P  P 

transformation  taking  P  to  the  origin  and  (A^>  =  [0,1]  *  [1,0]  .  Further 

let  $  =  J  (u  )  .  4>  will  be  called  a  standard  P-star  and  we  shall  call 

P  P  P  P 

members  of  S  =  { J  (A)!a€P  ,  AC“  }  the  standard  P  elements  of  4>  and 
P  P  P  P 

if  no  confusion  arises  denote  them  also  by  A  .  Note  that  4>  can  equal  4> 

P  Q 

for  P  4  Q  and  yet  5  ^  S  . 

P  Q 

LEMMA  4.3. 

i)  int  i>  is  a  domain. 

_ P _ _ _ _ 

ii)  There  are  not  more  than  Z  =  A(K)  possibilities  for  S  as  P  ranges 
-  p - 


over  R(V)  . 
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Proof . 


i)  Follows  Immediately  by  lemma  3.7, 


ii)  The  result  will  follow  if  we  can  show  that  there  is  a  finite  collection 

of  squares  S,,...,S  .  such  that  (J  (A)  |a€D  AC“  }  =  {S.....S  \  } 

1  -n(k)  p  1  p  1  ir(k) 

To  this  end  we  show  by  induction  that  the  vertices  of  each  J  (A)  must 

.  _  P 
_([/ a  yllnt+l) 

have  coordinates  of  the  form  (kn,^n)  where  n  =  2  v  82  J  k,£  integral 

|k[,j£|  £  2K^  and  that  diam  J^CA)  =  2^n  for  some  X  =  1,2,3,...  Here 

denotes  the  integral  part. 

Now  if  A  ,  ACtOp  we  can  construct  a  sequence 


A  =  A  , .  .  ,  A  =  A 
p  o  n 


having  the  properties  mentioned  in  Lemma  3.7. 

Clearly  any  vertex  of  J(A^)  satisfies  the  above  inductive  hypothesis. 

So  suppose  it  holds  for  J (A^) , . . . , J(A^)  0  £  i  n-1  .  Then  by  Lemma  3.1 

(3)  diam  =  2s  diam  A^  ,  s  integral,  giving  diam  =  2S  diam  J^CA^)  = 

2^+Sn  for  some  integral  A,s  .  Since  A . . , C  cu  we  have  as  in  4.2  diam  A...  > 

l+l  p  i+1  — 

diam  A  /k  .  This  gives  diam(J  (A...))  >  —  and  so  A  +  s  >  0  and  so  we 
p  p  l+l  —  K 

conclude  that  diam  J^(A^+^)  =  2%  with  t  >_  1  integral. 

Suppose  for  the  moment  that  A^+^  OL  diam  A^  .  Appealing  to  Lemma  3.1  (1) 
we  see  that  the  two  vertices  of  J  (ft,  which  are  the  end  points  of 

p  l+i 

-L(A  ,  i )  n  J  (A.)  must  satisfy  the  induction  hypotheses.  Since  J  (A,...) 

P  I'  1  p  1  p  It 1 

is  a  square,  it  follows  also  for  the  other  two  vertices.  The  bound  on 
|k|,|£|  is  a  consequence  of  diam  <t  £  2K  .  The  case  diam  A^  <_  diam  A^+j 
follows  by  a  similar  argument. 

COROLLARY:  There  is  not  more  than  Z(K)  of  different  possible  domains  $ 


w 
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Given  &€V  we  denote 

(4.3)  Q* (A)  =  Uk>  )p€R(P)  ,  ACa)p} 

As  in  the  proof  of  Lemma  4.1  we  have  Q*(A)£Q(A),  where  Q(A)  is  a  square 
with  the  center  at  the  middle  of  A  and 

(4.4)  diam  Q(A)  2(K-l/2)  diam  A 

LEMMA  4.4 

i)  If  A',A"^P  and  A',A"CQ*(A)  then 


_1 

(4.5)  _ diam  A  '  >  dram  A 

and 

(4.6)  diam  A'  -kf  diam  A" 

_ _ _  _  K" _ 

ii)  Int  Q*(A)  is  a  domain  . 

iii)  If  P€R(V)  and  wpCQ*(A)  then  on  A 

(4.7)  j  D^v  j  <  K(diam  A)  ^  . 

P  ~ 
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Proof . 


hence 


Let  A ' C  Q* ( A) 


then  for  some  P'  €R(P) 


we  have 


A,A'  €  o> 


(4.9) 


diam  A' 


>  —  dlam  w  , 
—  K  p' 


^  —  diam  A 


Further  we  have  for  A"coj  „ 

P 


(4.10)  diam  A  >  diam  w  ,,  >  diam  A" 
—  K  p  —  K. 


(4.9)  and  combination  of  (4.9)  and  (4.10)  yield  (4.5)  and  (4.6). 

ii)  We  have  to  prove  only  that  Q*(A)  is  connected.  This  follows 
immediately  from  the  Lemma  3.7. 


iii)  From  Lemma  3.5  (ii)  and  (iii)  we  have  0  <_  v  <_  1  and  so  obviously 


and 


!l)1v 


[  min  diam  A]  ^ 

AC  oj 
P 


[  min  diam  A 

AC  to 
P 


n 


(4.9) yields  the  Lemma. 

so 

Let  K  <  2  .  Denote 


Q**(A) 


(x  € Q*(A) |dist(x,3Q*-3fi) 


diam  A 
-  Tq+T 
2  ° 


} 


LEMMA  4.5. 


i)  Let  'JCQ*(A)  .  Then  P€Q**(A)  . 
_ _ _ P_  _ _ 

li)  ACQ**  (A) 
iii)  int  Q**(A)  is  a  domain. 


Proof . 


s„+2 


i)  Assume  that  P^Q**(A)  then  dist  (P,3Q*-3£i)  <  diam  A/2  ° 

Because  v^  =  0  on  3Q*-3H  ,  (4.7)  leads  obviously  to  a  contradiction. 

ii)  Assume  that  A<^Q**(A)  .  Then  there  is  a  vertex  P*  of  A  such 

that  dist  (P* ,  3Q*-3fi)  <  (diam  A)2  ^s°+^  .  Obviously  w  =  £  v  =  1 

a>jPQ*(0)  P 

on  A  .  By  the  same  argument  as  leading  to  (4.7)  we  see  that  |D^w|  <_  K(diam  A) 

Because  w  =  0  on  3Q*-3ft  we  have  a  contradiction. 

iii)  Let  x£Q**(A)  ,  x€A  C  oj  .  Then  by  Lemma  3.7  there  exists  sequence 

o  p 

n 

A  A  =  A  such  that  V  =int  U  A.  is  a  domain.  Because  diam  A.  > 

°  1  i-0  1  J  ■ 

—  diam  A  by  (4.2)  it  is  readily  seen  that  int  Q**flV  is  a  domain.  This  leads 
K 

immediately  to  the  desired  result. 

Let  if*  =  Jp(Q*(A))  (analogously  as  <fp)  and  if**  =  Jp(Q**(A)).  Then  we  have 


LEMMA  4.6.  There  exists  not  more  than  Z*(k)  (respectively  Z**(k)) 
possible  domains  int  if*  (respectively  int  if**)  . 


Proof.  The  first  part  of  the  lemma  follows  from  Lemma  4.3  and  is  a 
corailary.  The  second  part  follows  by  analogous  argument  when  realizing  that 
if**  is  composed  by  squares  with  diam  n/4  when  n  was  introduced  in  the  proof 


of  Lemma  4.3. 
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5.  THE  APPROXIMATION  PROPERTIES  OF  M(P) 


Ln  this  section  we  will  analyze  the  approximation  properties  of  M(D)  . 
Let  P.  be  the  mapping  of  C°(0)  onto  M(D)  such  that  (IIu)  (P)  -  u(P)  , 
p€R(D)  .  By  Lemma  3.5  we  can  write 

n  u  =  y  u(p)v 

P€R(V)  P 

Further  define  the  operator  ^  ^  mapping  C  (fi)  into  M(P)  by 

(5.1)  nQ*(A)(u)  =  )'  u(P)v 

P€R(D)  P 

io  CQ*( '.) 

Clearly  supp ( 1'^  ^u^(u))CQ*(A)  and 

(5.2)  r!Q*(A)  (u)  =  I'  (u)  on  A  . 


For  given  A€0  we  define  as  the  invertible  affine  transformation 

2  2 

of  TR  onto  IR  taking  A  into  standard  unit  square  [0,1]  *  [0,1]  =  S  . 

LEMMA  5.1.  There  exists  a  constant  C  dependent  only  on  K  such  that 
for  any  w€H^(lQ  and  £=0,1 


(5.3) 


I  wo  J, 


(nQ*(AW 


i  p 

H4(S) 


<  C I  I woj 


-li 


H^(JaQ**(A)) 
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Proof.  By  Sobolev  imbedding  theorem  we  have 


(5.4) 


|  woJ 


-1, 


c°(Jaq**(a>) 


<  C I  I woJ 


-1, 


H2(JaQ**(A)) 


Because  of  Lemma  4.6  there  are  at  most  ?.**(K)  different  domains  J^Q**(A)  , 
where  Z**(K)  depends  on  K  only.  This  shows  that  C  in  (5.4)  depends  on 
K  only. 

Using  (5.1),  Lemma  4.4  and  Lemma  4.5  i)  we  get 


(5.5) 


(n 


Q*(A) 


w)°j" 


I  l 
H  vS ) 


<  C I  I woJ 


-li 


H2(JAQ**(A)) 


where  C  depends  on  K  only. 

Since  by  Lemma  4.5  ACQ**(A) 


(5.6) 


H^(s) 


< 


I  IwoJ-1! | 

H^(JaQ**(A)) 


Combining  (5.5)  and  (5.6)  we  get  the  Lemma. 


LEMMA  5.2.  There  exists  a  constant  C  (dependent  only  on  K)  such  that 


for  any  w€H  (fi)  and  £  =  0,1 


(5.7)  IwoJ"1  -  (IIQ*(A)w)oJ  X|  iClwor1) 

H^(s)  A  HZ(JAQ**(A)) 

Proof.  Suppose  z  is  any  function  bilinear  on  Q**  .  Then  on  A  we 
Q*(A) 

have  z  =  z  and  therefore 
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(w+z)  «J  1  -  (!IQ*(A)  (w+z)) °J  1  =  woj  1 
A  A  A 


oi^W1 

A 


and  hence  using  lemira  5.1  we  get 


Iwoj"1  -  ( nQ*(A w)°J  1 


■1  .i 


A  m:  w/  -'a  l  1  Q  --  ^  lW°^A  “Zl  I  ? 

A  Hy'(s)  z  A  H (J  Q** (A) ) 

A 


where  z  is  an  arbitrary  bilinear  function  on  J^Q**(A) 


Because  by  lemma  4.6  there  is  only  a  finite  number  (depending  on  K  only) 
of  different  J^Q**(A)  ,  and  by  lemma  4.5  any  int  J^Q**(A)  is  a  domain,  we 
have  (see  e.g.  (  2  1,  pp.  )  that 


Inf[(w°J  *-z[(  ?  <  c|w<J 

ft(Jaq**(a))  h/(Jaq**(a)) 

The  theorem  follows  immediately. 

Now  we  have 

THEOREM  5.3.  There  exists  a  constant  C  depending  only  on  K  such  that 
for  any  A  CP,  u€H^(li)  and  l  =  0,1 

(5.8)  I  1 11-rr  u  I  I  <  C(diam  a)2  ?  |  |  u  |  | 

H  ' (A)  H  (Q**(A)) 


Proof .  Using  (5.2)  we  obtain  (5.8)  from  lemma  5.2  by  standard  scaling 
argument . 


THEOREM  5.4.  Let  uCH2(fi)  then  for  l  =  0,1 


with  C  depending  only  on  K  . 


Proof.  We  have 

I |u-nu| |2  =  I  i!u-nu|!2  £  c  H  |u||2  h2(2_fc)(P)< 

H  (ft)  A€0  HX(A)  A€fl  H  tQ**(A) ) 

1  C  l  ||u||2  h2(2"^(P)  <  C  l  A(A)|  |u|  |2  h2(2_i°(0) 

b€V  r(Q*(A))  h€V  HZ(A) 

where  A(A)  is  the  number  of  Q*( A’)  such  that  ACQ*(A')  .  Using  lemma  4.1 
we  see  that  X(A)  £  MN  £  C(k)  and  so  theorem  5.4  is  proven. 

2 

Remark:  In  Theorems  5.3  and  5.4,  the  restriction  u^H  (ft)  can  be  weakened 
to  u€H2(A)  for  every  A€P  . 

Theorem  5.4  shows  that  M(V)  has  the  same  basic  "interpolation"  properties 
as  the  usual  finite  element  spaces.  The  spaces  M(P)  are  more  flexible  than  the 
usual  spaces  defined  on  quadrilateral  meshes.  The  space  M(D)  allows  us  to 
make  a  refinement  and  still  keep  square  elements.  The  restriction  to  K-meshes 
is  from  a  practical  point  not  important.  A  more  essential  restriction  is  that 
we  deal  only  with  squares.  How  to  overcome  this  resctiction  with  respect  to  the 
implementation  and  the  theory  will  be  discussed  elsewhere. 

The  use  of  the  spaces  M(V)  does  seem  to  have  a  major  advantage  over  tri¬ 
angular  elements  because  of  programming  and  data  management  simplicity,  expecially 
when  some  form  of  automatic  or  adaptive  mesh  generator  is  envisaged.  One  mani¬ 
festation  of  this  is  that  the  element  can  be  uniquely  defined  by  the  binary 
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expansion  of  the  coordinates  of  its  center,  the  length  of  such  an  expansion 
indicating  the  size  of  the  element.  Since  all  the  elements  have  moduli,  a 
scaling  factor  the  same  geometric  shape,  the  calculation  of  the  stiffness 
matrix,  etc.  is  simplified.  In  addition  this  seems  to  be  important  for  the 
practical  effectiveness  of  the  error  estimation. 

We  shall  analyze  now  the  approximation  properties  of  M(P)  when  u€H'*'(51). 

2 

Before  being  that,  we  introduce  some  notations.  By  p(x),  x ClR  we  denote 
a  molifier,  a  function  with  all  derivatives  continuous  p(x)  >_  0  ,  p(x)  =  0 
for  |xl  >  1  ,  u(0)  =  1  ,  and 


p(x)dx  =  A  .  For  e  >  0  ,  let  y^(x) 


IR„ 


f.  2A 


h(x/e) 


Let  JICIR^  be  an  arbitrary  bounded  domain  and  RP  its  p-neighborhood . 


For  u  defined  on  R  ,  we  put 


u  =  u*p  ,  c  <  p 

c  e 


Obviously  u  is  defined  on  R  . 

J  e 

2  , 

Further  for  any  t  6TR  ,  |t|  <_  1,  e  <  p  let 


Ute(x)  =  u(x+tt) 


Then  utC  is  also  defined  on  R  ,  and  we  have 


(5.10)  1  I  u 


E  '  V(R) 


lCC_1||u|i 

H  (RP) 


lUf-Ul  I  a  1  Cf;1  ?'l  lul  I  1  _  »  O'  *  0,1 

H  (R)  H  (RP) 
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|utE-u| |  1  Ce1_2| |u|| 

H  (fi)  H  (ST) 


with  C  independent  on  ft  ,  u  and  e 


Let  us  prove  now 

LEMMA  5.5.  For  every  P€R(V)  ,  let  a  function  w  6H2(ft)  be  given.  Let 
-  -  -  -  _ _ _ E - 

A  €  V  and  A C  w  .  If 

Po 


w  =  n  7  v  w 

P€R(P)  p  p 


then  for  l  =  0,1 


(5.11) 


h£(a) 


C[(diam  A) 2  ^ 


wp  I  I  2 

Po  H  (Q**(A)> 


+  /  y  |  |--w  +w  |  |  .  (diam  A)  * 

P€R(P)i=0  Po  P  H  1(Q**(A) ) 

A€u 

P 


Proof.  On  A  we  have 

w  =  n  y  vw  +  n  y  v[-w+w] 

P€R(5)  P  Po  P€R(A)  p  Po  P 

=  nQ*(A)w  +  nQ*(A'  l  v  ( -w  +w  ] 
po  P€Ri A)  p  Po  P 


Using  lemma  4.4  and  Leibnitz's  rule,  we  have  for  l  =  0,1,2 
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(3.12) 


l  I 

A '  n  Q**(A)f*0 


(q**(a)  n  a'  ) 


C  )  | |-w  +w  | !  .  (diam  A) 

i=0  po  P 


i-P 


Therefore  on  A 


„  _  nQ*(A) 
-w  =  w  -li  w 


/  V  [-W  +w  ] 

P€R(A)  P  Po  P 


[nQ^A)( 


y  v  f-w  +w  1 
P€R(A)  p  po  P 


)  V  [  -W  +W  I  ) 

P€R(A)  P  Po  P 


Applying  ttieorem  5.3  and  noting  the  remark  after  theorem  5.4,  we  get  the  lemma. 

LEMMA  5.6.  Let  v€H1(IR2).  Associate  to  every  P€R(V)  a  sector 
t  €  IR  ,  | t  I  r  1  and  a  number  X  such  that 

_  P  -  - _ _ 


>, 

P 


mi n(diam  A  )' 

A€  , 

P 


1 

K.512 


Let  further 


t  3 

,  P  P, 


1 

-  i> 
64 


and 


w  =  n  y  v  w 

P€R(V)  p  p 
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Then  for  any  A€P  and  i,  =  0,1 


|v-W| I  £  C(diam  A)1  l\ |v| | 

H  (A)  H  (Q**(A))  P 


(5.13)  p  =  max  A 


mpCQ*(A) 


Proof.  Let  Adw  and  ACw  .  Then  using  lemma  4.4  we  have 

-  pp  ° 

ro 

Ap,Ap  £  C(K)  diam  A.  Using  (5.10)  we  get  for  i  =  0,1,2 


w  -w  . 
P  p  1  '  i 


<  C(diam  A) 1  i I | v| 


o  K  (Q**(A) ) 


H1([Q**(A)]2p) 


By  lemma  5.5  we  get  readily  the  lemma. 


LEMMA  5.7.  Let  p  be  defined  by  (5.13),  then 


[Q**(A)]4d  DnCQ*(A) 


The  lemma  follows  easily  from  the  definition. 

THEOREM  5.8.  Let  u€H^(fl)  respectively  H^,r(fl)  .  Then  for  any  V 

1  r 

with  h(0)  _<  there  exists  w€M(V)  respectively  M(P)f"lH  ’  (ft)  such  that 

for  all  A€P  and  i  =  0,1 

(5.14)  [  | ii— w |  |  <  C(diam  A)1  ^  |  |u|  |  , 

H*(A)  H  (Q*(A)) 


1 
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with  C  independent  of  V,  u  and  A  . 


Proof .  The  extension  of  u  onto  a  neighborhood  ft0  ,  has  the  properties 

2po 

listed  in  theorem  2.1.  Denote  by  [Q**(A)]  the  2Po-neighborhood  of  Q**(A) 

p  =  min[p,p/2Sl  where  p  is  defined  bv  (5.13)  and  8  was  introduced  in 
o 

Theorem  2.1.  Then  bv  lemma  5.7  and  theorem  2.1 


IN!  x  2p  -  I  i  4p  -  cl  M  i  i 

ii  ( [Q**(A)  j  )  Hi<[Q**(A)]  PDft)  Ha[Q*(A)] 

Let  u€H^(ft).  Select  t  =  0  in  Lemma  5.6  for  all  P€R(V)  .  Then 

P 

for  W  constructed  in  Lemma  5.6  we  have  W€M(0)  and 


|u-W| |  <  C| |u| 

Hl(A) 


!  !  2p  i  c  I  1 u  |  |  , 


and  the  first  part  of  the  theorem  is  proven. 

1  F 

We  have  now  to  show  that  when  u€H  ’  (ft)  then  we  can  choose  W  so 

that  W(P)  =  0  for  all  P€R(V),  P€F  .  To  every  P€R(P),  P^F  ,  we  take 

t  =  0  .  Tf  P  €  3ft ,  P  €  F  ,  p€?V  then  we  take  t  the  outward  unit  normal 
P  P 

(if  P  is  a  corner  of  ft  then  the  normal  is  bisecting  the  outside  angle). 

If  P€lf’  then  it  is  easy  to  see  that  we  can  select  vector  t  (pointing  in 

to  the  sector  ftPy.  in  Fig.  2.1)  and  not  necessarily  of  the  unit  length  such 

that  w  (P)  =  0  also.  This  finishes  the  proof. 

P 

Finally  we  prove 

i  °  I  i  r 

THEOREM  5.9.  Let  u€H  (ft)  [respectively  H  (ft),  H  ’  (ft))].  Then  there 


exists  w€M(V)  respectively  M (17)  fl  H 3  (S2)  respectively  M(P)flH1,r(ft)  such 


that 
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(5.15)  l  j  j v  (u-w)  |  | 2  <c||uj|2 

P€R(V)  P  H  (SI)  H  (fi) 


with  the  constant  C  dependent  only  on  K  . 


Proof .  By  Theorem  5.8  we  see  that  for  any  A€V  and  1  =  0,1 


(5.16)  ||u-w||  £  C(diam)1  ^||u|| 

H  (A)  H  (Q*(A) ) 


and  applying  Leibnitz's  rule  and  lemma  4.4,  we  obtain 


(5.17)  | | v  (u-w) | |  <C||u||  . 

P  Hi(A)  n  (Q* (A) ) 

for  all  P€R(V )  such  that  w^CQ*(A)  .  By  lemma  4.1  there  are  not  more  than 
M(K)  such  nodes  so  we  have 

l  I |v  (u-w) I |2X  <  l 

P€R(V)  P  H  (A)  P€R(P) 

w  CQ*(A) 


and  therefore 


(5.18)  l  |  |  v  (u-w)||2  <  C(K)  l  IjuH2. 

P€R(P)  p  H  (fi)  ts€D  H  (Q*(A)) 

By  lemma  4.1  there  are  not  more  than  C(K)  different  A'  €.V  such  that 
Q*(A’)3A  and  so  (5.18)  yields  (5.15). 


v  (u-w)  |  |2  <_  C(K)  |  |  u  |  |  2 

P  H  (A)  n(Q*  (A) ) 
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Remark  to  Theorem  5.9.  Assume  now  that  int  supp  u  =  Q*  .  Then 


supp  w  =  U  w 

P€R(V)  P 

o-  fUl*^0 
P 


This  observation  follows  immediately  from  (5.16). 


Remark:  Theorem  5.9  is  closely  related  to  (3.1)  of  [  9  ]  which  is  an 

essential  part  of  the  a-posteriori  error  analysis. 
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6.  THE  MODEL  PROBLEM,  ITS  FINITE  ELEMENT  SOLUTION 


AND  THE  BASIC  A-POSTERIORI  ESTIMATE 


As  a  model  problem  we  will  discuss  the  case  of  plane  elasticity  for  a 

body,  homogeneous  and  isotropic  on  every  Q  making  up  the  domain  ft  . 

0 

Let  H1(ft)CH.CH1(ft)  i  =  1,2  where  H.  =  H1,ri(ft)  as  defined  in 
Section  2. 

Let  Hq  =  H^  x  H£  ,  with  u  =  (u^,u2)  and  consider  on  Hq  x  Hq  the 
following  bilinear  form 


(6.1)  B(u,v)  =  B(u1,u9;v1,v2)  = 


ft 


r  3u  3v  3u„  3v 

[(A+2^ (  3^  3^  +  ir2 1Z?  + 


3u  3u„  3v  3v  3u1  3v  3u  3v 


5-} 


We  assume  that  A  and  p  are  positive  and  constant  on  every  Q  of  ft  . 
the  constants  A  ,  u  are  the  usual  Lame  constants. 

We  will  assume  that  there  exist  constants  0  <  crc2  <  00  such  that 
for  any  u  €  H^ 


(6.2)  C  |  | u |  |  1  B(u,u)  £  C J  |u|  | 

H 1  (ft)  Z  Hl(ft) 


wi  th 


2 

HX(n) 


H1  (ft) 


Ik 


2 

H1(n) 


V 
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Obvious  Iv  B  ( u ,  v )  ‘  =  B(v,u)  and  |B(ti,v)|  <  C.||u||  ||v||  •  There- 

H  (SO  H  (ft) 

fore  on  H^,  B(u,v)  is  a  scalar  product  with  the  energy  norm 

(6.3)  j  | |uj  j  j2  =  B(u,u) 

(6.2)  shows  that  the  energv  norm  |  |  !u|  I  I*'  Is  equivalent  with  |  j  u  j  j  . 

H  <fi) 

The  problem  P(H  w,g),  w  C 11 1  (S'O  x  H1  (U)  ,  g€ll°(fl)  x  H°(S2)  consists 
of  finding  x  H  ^  ( ..2 )  such  that  u  -  w€Hq  and 

(6.4)  B  (  u  ,  v )  =  (g.v)  ,  WvCH^ 

where  we  have  written 


(g,v) 


(fi]v1  +  g2v2)dx 


It  follows  by  the  standard  theory  that  u  exists  and  is  uniquely  determined. 

The  function  u  wiLl  he  called  the  exact  solution  of  the  problem  P  . 

Let  now  M.  (P)  =  .U(P)  0  H  (::),  i  -  1,2  and  M0(D)  =  M^P)  x  M2(D)  . 

■) 

Assuming  that  wC  (M(P))‘L  in  the  problem  P(H  ,w,g)  ,  the  finite  element 
solution  H  of  P  with  respect  to  M  consists  of  finding  U€M(P)  x  M(P) 
such  that  U-w^M  (P)  and 


(6.5)  B (U , v)  *  (g.v),  Vv€M0(P) 


Just  as  for  the  exact  solution  it  follows  that  the  finite  element  solution  U 
exists  and  is  uniquely  determined. 


39 


Finally  we  denote  by  e  =  U  -  u  the  error  of  the  finite  element  solution 


We  will  be  interested  in  an  a-posteriori  estimate  of  the  norm 


| '  jej  | |  or  some  norm  equivalent  to  it.  We  will  design  an  estimator  E 
depending  only  on  the  known  finite  element  solution  —  which  will  be  related 
to  the  error  norm.  As  estimator  E  is  called  an  upper  respective  to  lower 
estimator  if  there  exist  constants  respectively  A^  independent  of  V 

and  u  ,  U  such  that 


IMI  i  AyE 

respective 

AlE  <  |le|| 


THEOREM  6.1.  Let  wCM(P)  u  be  the  exact  solution  of  the  problem 


P(Ho>w,g)  and  U  its  finite  element  solution  with  respect  to  M  .  Then 
there  exist  strictly  positive  constants  C^,  depending  only  on  K  , 

u  ,  A  and  0  such  that 


(6.6) 


"0 


JL  B(nn*nn)  I 


PCR(O) 


<  C1  I  Mn  ) 

“  1  P€R(V)  p  p 


where 


(6.7) 


n  € H ( ii'  )  =  {v€Hnjv=0  on  fl-u>  } 
P  P  O'  p 


and 


(6.8)  B(n  ,v)  =  B(e,v) ,  \fv€H(w  ) 
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Proof.  It  follows  from  the  definition  of  U  and  u  that 


B  ( o  ;  v )  =  0,  Vvfjll 


and  therefore 


B(e,e)  =  B(e ,e-v)  =  B(e,  }  v  (e-v))  =  £  B(n  ,v  (e-v)) 

P€R(P)  p  V€RTV)  p  p 


Using  Schwarz's  inequality  we  get 


B  (  e  ,  e ) 


P€7 1(0) 


in  i  | 

P 


!  I 


v  (e-v)  !  |  I 
P 


t  1 

?€R(V) 


2 


1/2 

] 


By  (6.2) 


P€R(V) 


v  (e-v 
P 


(e-v)  |||“  C|  |  v  (e-v)  |  !“1 


H  (Si) 


and  therefore  by  Theorem  S.9  for  a  proper  choice  of  v 


L 

P  €i(V) 


i  v  ( e-  v ) 
P 


i  2 


Cl  iel 


H1  C«) 


2 


Hence  we  have 


P€R(P) 


wh  i  cii  proves  right  hand  side  of  (6.6). 
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Let  us  prove  now  the  left  hand  side  of  (6.6).  Define 


w . 
J 


j  =  1 


where 


v.  are  the  sets  introduced  in  lemma  4.2.  We  have  now 
J 

(6.9)  B ( e , w . )  =  l  B(e,n  )  =  j  B(n  ,n) 

J  P£,  P  pgx  P  P 

j  'j 

Because  for  P  ^  Q  ,  P,Q£x.  jnt  (supp  n  )  Hint  (supp  rig)  =  0  we  have 

B(n  ,  n,J  =  0  3nd  hence 
P  Q 

(6.10)  B(v.,w.)  =  l  B(n  ,n  ) 

■>  J  I  ‘ 


Further 


|B(c. 


i  i 


and  hence  using  (6.9)  and  (6.10)  we  get 


I  I  l”.l  I  !Z  -  y  III"  1112  <  IIHII  III-, III 

P€X.  P  J 


and  hence 


|W 


P€X: 


Because 


j  ranges  over  1  I  we  get 


2 


(6. 11) 


>:  Mini 

P  €R(V)  p 


By 


lemma  4.2  -I  <  L(k)  and  therefore 
We  have  proven  in  Theorem  6  that 


(6.11)  proves  the  left  side  of 
the  estimator 


(6.6). 


t  (U)  =  l  B(n  , 

P€R(V )  p 


nP) 


is  simultaneously  on  upper  and  lower  estimators.  The  individual  terms  Hp 
are  determined  locally  on  separate  stars  w  .  Let  us  underline  that  al¬ 
though  the  unknown  error  e  is  present  in  the  definition  of  np  in  (6.8) 
we  have  not  to  know  it.  In  fact 


B(e,v)  =  B(U-u,v)  =  B(U,v)  -  B(u,v)  =  B(U,v)  -  (q,v) 

Remark.  The  proof  of  the  theorem  6.1  follows  very  closely  the  ideas 
in  [  9  ]  . 

We  assumed  in  Theorem  6.1  that  w€M(V)  .  It  is  obvious  that 
w  =  (wpiv^)  influences  the  solution  only  by  its  values  at  the  boundary 

3L!  ,  more  exactly  on  F  .  In  general  when  w^M(P)  we  replace  w  by 

<\,  __ 

wCM(P)  and  estimate  the  norm  of  the  solution  of  the  problem  P(Ho,w-w,0)  . 

Usually  it  is  easy  to  explicitly  construct  a  function  z£  H^(ft)  x  H^(ft) 

=  (w-w)  j  on  ,  i  =  1,2  and  the  desired  estimate  is  then  simply  |||z||| 

In  practical  cases  we  can  expect  that  lllz||  I  is  much  smaller  than  |||e|||  . 
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7.  THE  A-P0STERI0R1  ERROR  ESTIMATE 


Let  us  denote  by  Q  any  of  the  squares  comprising  the  domain  £1  and 
let  Oq  be  the  unit  square  [0,1]  x  [0,1]  .  Assume  that  €H^(Qo>  , 
i  =  l,...,n  are  given  and  denote  by  Z  the  linear  span  of  ,  i  =  l,...n. 

Definition  7.1  Let  p  >  0  ,  e  >  0  .  By  ^(Z.p.e.Q)  we  denote  the 
family  of  functions  t,€  H°(Q)  such  that  for  any  square  S  =  [a^,a^+h]  x 
[a2»a2+h]CQ  the  following  properties  are  fulfilled 

a)  C€H1(S) 


E)  There  exists  =  z(h  ^(x-a)),  z€Z  and  constant  M  (both  depending 


on  C  )  such  that 


i)  I  U-S  !  I  I  M-h 

H  (S) 


ii)  U-U  ,  £  m 

HL(S) 


iii)  U  |  ,  >  pMh 

H  (S) 


Let  us  illustrate  our  definition  by  a  few  examples. 

(i)  Let  Z  be  the  set  of  all  polynomials  of  degree  less  than  or  equal 

to  n  .  Then  any  polynomial  of  degree  £  n  on  Q  belongs  to  the  family 
H(Z,p,h,Q)  with  p  and  r  arbitrary. 

ii)  Let  Z  be  the  set  of  guadratic  polynomials.  Then  £  =  sin  x^ 

belongs  to  the  family  H(Z,p,e,Q)  for  e  £  1  and  some  suitably  chosen  P  . 
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iii)  The  family  H  is  typically  characterized  by  2  being  the  poly¬ 
nomials  of  degree  m  ,  say,  and  then  we  take  to  be  a  suitable  Taylor 

expansion  of  f  ,  Then  (i)  and  (ii)  are  more  or  less  standard,  and  (iii) 
states  that  is  not  'degenerate". 

LEMMA  7 . 1 .  Suppose  that  f (2 ,p ,g ,Q)  then  there  exists  C  (depen- 
den  t  on  (Z,p,e,Q),  such  that  for  any  square  SCQ 

(7.1)  | f  |  <  C(diam  S)"1  inf  ||f-d|| 

H  (S)  d  =  constant  H  (S) 

functions  on  S 


Proof .  1)  Define  F€H^(Q^)  by 


F(x)  =  f(a  +  hx) 

Now  by  the  assumption  there  exists  F  (x)€7.  such  that 

o 

(7.2)  j  | F-F  ( |  <  M  , 

H°(Qq) 

(7.3)  | F-F  |  <M  , 

°  “  <V 

(7.4)  If  I  ,  >  pMh_E  . 

0  H1(Q  )  " 

o 


Denote  now  by  F  respectively  F^  the  average  of  F  respective  to  Fq  on 

Q  .  Then  ve  have 
o 


(F-F  )-(F-F  ) |1  <  || F-F  | |  <  M 

H°(Q  )  °  H°(Q  )  “ 

o  o 


(7.5) 
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because  F-F0  is  t,le  H°  projection  onto  the  set  of  constant  functions  on 

%  • 

Because  Z  is  finite  dimensional  space  there  exists  C  (dependent  on 
Z  )  such  that 


(76)  If  I  =  If  -f  I  ,  <  c | | f  -f  | | 

0  H1(Qo)  °  °  H1(Qq)  °  °  H°(Qo) 


and  hence  combining  (7-4),  (7.5)  and  (7.6)  we  get 


M  < 


,  £  “I  I 
h  p 


°  H1(Qq) 


<  Ch 


:P_1| IF  -F 


H°(Qo) 


and 


•  I CF-Fo)-(F-Fo) 


H°(Qo) 


Cp~1h€ 


|  |f  -f  | 

1  '  o  o  1 


H°(Qo) 


On  the  other  hand 


F_F  H°(Q  )  "  F°"F°  H°(Q  )  "  HW(QJ 

o  o 


||(F-Fo)-(F-Fo)|| 


and  hence  for  Cp  ^h£  <  1/2  we  have 


(7.7)  I  I F-F | | 


>  i  F  -F 

H°(Qo)  “2  °  °  B°(Q0) 


Further  from  (7.3)  and  (7.4) 
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(7.8)  | F I  <  I F  |  +  |F  -F|  .  c  | F  |  . 

HX(Q  )  °  H 1 (Q  )  °  H1(Q  )  °  H2(Q  ) 

o  o  o  o 


+  M  <  IF  I  (1  +  p'V)  <  2 1 F  I 

Hi(Q  )  0  H1(Q  ) 

o  o 


for  hl'p  ^  <  1 


Therefore  by  (7.6),  (7.7)  and  (7.8),  we  have 


IFI  .  <  2|Fo|  <  2c||Fo-FoM  o  <  4C  |  |  F-F  | 


H1(Q  ) 
o 


H  <v 


h  <V 


H° (Qq) 


so  upon  rescaling  back  to  S  (7.1)  is  proven  in  the  case  diam  S  £  ho(p,e,Z) 
2)  Suppose  now  that  diam  S  >  h  (p,e,Z).  Put 


r,.  „-iINT 
[diam  S  .  , 

a  ■  K~J  + 1 


1  NT 

where  [']  denotes  as  before  the  integral  part.  Clearly  we  can  divide 

2  2 

S  into  o  congruent  squares  S^,  i  =  l,...,o  with 


, .  c.  diam  S  , 

diam  S.  =  — • — —  <  h 
l  ct  —  o 


Thus 


!<  =  l 

H  (S)  1=1 


2 

,|2  <c(^fLS)  l  inf  |  |f-d  i  |2 

H  (S  )  i=l  H°(S.) 

o  , _ .  .  i 


d=constant 
function 


<  Co2 (diam  S)  2 


inf  | | f-d | 
d=constant 
function 


H°(S) 
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But  o  <  1  +  — - ^  <  C  and  the  result  follows. 

—  n  ^ 

o 

k  k 

We  introduce  the  family  H(Z,p,e#Q)  for  Q  =  with  0~  as  in  section 

b  b 

k 

2.  Because  fi  consists  of  a  finite  number  of  Q.  it  is  clear  that  we  can 

b 

k 

extend  the  family  H  into  H(Z,p,e,ft)  so  that  the  restriction  on  Q  =  Q_ 

b 

is  the  family  fl(Z,p,e,fi)  . 

In  what  follows  we  will  assume  that  we  are  concerned  with  problems  P(Ho>w,g) 

introduced  in  section  6  where  g  H  (g^,g2)  ,  g^  € H (Z ,p ,e ,0)  . 

We  shall  discuss  now  the  error  estimate  of  section  6  in  more  detail.  We 

have  shown  in  theorem  6.1  that  the  essential  part  of  the  estimate  is  the  norm 

of  n  which  is  defined  on  the  star  u>  .  In  section  4  we  introduced  the 
P  P 

standard  star  4>  (see  lemma  4.3)  as  the  image  of  to  under  the  mapping 

J  . 

P 

On  we  will  define  now  the  space  H  =  H,  x  H„  with 

P  P  l.P  2,p 


H  (*  )  -  tv  -  (v,  ,v?)  €H1(4>n)  x  H1(4’n)|voJ  €H(w  )} 

PP  ±  £.  P  PPP 


^tb  fif(w  )  defined  in  (6.7). 

In  6.1  we  defined  the  bilinear  form  B 


Let  us  define  the  form  B 


defined  on  x  with  the  same  expression  as  in  (6.1)  but  with  integra- 

*  -  _  ^  i 

tion  over  ,  and  X,  u  instead  of  X,  y,  X  =  X<>J  ,  m  =  U°J 


LEMMA  7.2.  The  bilinear  form  B  is  such  that  for  any  u€H  , 
- -  •  ___  _ _ P _ _ _ P 


(7.9)  C  |  |u||2  <  B  (u,u)  lC2||u||2 

1  HL($  )  P  H*(*  ) 


with  constants  ,  i  =  1,2  dependent  only  on  X,y  and  K  but  independent 
of  P  . 


48 


Proof.  For  each  P€R(P)  there  are  line  segments  T'’  C.  3<{>^  (i  =  1,2) 

both  being  edges  of  standard  P  elements  of  <S>p  and  such  that  u  €  Hp  =  0 

on  T(i)  .  Applying  now  lemma  4.3  we  see  that  there  are  only  finite  number 

of  different  cases  of  the  domains  <J>  and  the  line  segments  F^^  for  each 

P 

the  Korn  inequality  holds.  Therefore  (7.9)  holds  with  |-|  .  replacing 

"  <V 

| | ■ | |  .  .  Because  these  two  norms  are  equivalent  with  constants  depending 

H  <*p) 

on  (l>p,r(1)  .r^2-*)  we  get  (7.9)  immediately  when  using  again  lemma  4.3. 

Denote  now  by  M  (f>  )  the  set  of  all  u€H  ($  )  being  bilinear  on 
P  P  P  P 

every  standard  P  element  of  the  standard  star  $  .  Further  we  denote 

P 


-  _1  I  2 

e  =  e°J  ,  u  =  u°J  ,  U  =  U»J  ,  gi  =  g^J  ,  g£  =  gi(diam  A  ) 


i  =  1,2 


where  e,  u,  U,  g  were  defined  in  section  6. 


We  have 


B  (e  v)  =  0  ,  Vv  (<J>  ) 
P  P  P 


8p(e,v)  =  Bp(U,v)  -  Bp(u,v) 


|(“'V)  =  l  ^i.v>Lo ($  y  Vv€Hr 
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By  integration  by  parts  we  get 


Bjl’.v)  =  l  l 


n 


A£5  i=l  j  =  l 
P 


o .  ,  (U)n.v.ds 

1 ,  j  1  j 


.i(U)v.dx j 


3A  A 

where  (n^,n2)  is  the  unit  normal  on  each  edge  of  3A  pointing  outward  from 
A  and 


„  9U  3U 

°i,j  =  X(9^  +  3^}  +P 

„  ,  „  ,  32u„ 


f  aO  au.  1 

S7+a^J 


Li<U)  -  (»+„)  3^ 

.  .  .  .  22U 

4<u>  ■  n^sr2 


Let  F  (respective  Y0  )  be  the  union  of  all  vertical  ^horizontal) 

1  ,  P  ^  9  P 


edges  of  all  A  € ,  Then 


(7.10)  B  (e,v)  =  l  f-(fg  +L  (U)],v  ) 

P  i~l  L 


W  +  j-i  U‘<<’i.J<’,)>lVj).L2<riip)] 


where  indicates  the  jump  in  a  function  across  .  The  jump  has 

an  obvious  sense  if  the  relevant  edge  T  of  A  is  inside  4>  .  If 

P 

rC  3$  we  have  to  distinguish  the  case  of  whether  v.€H,  =>v.  =  0  on  T 

P  J  1»P  J 

or  not.  If  all  v^  =  0  on  T  we  will  take  the  jump  to  be  zero.  If  it  is 

not  the  case,  then  is  set  equal  to  zero  outside  ♦  ,  and  the  jump 

is  taken  the  usual  sense.  This  convention  will  be  used  through  the  paper. 

Because  U  is  bilinear  on  every  A  6 $  and  A,p  are  constant  on  A 

P 


we  see 


that  L^(U)  are  constant  on  every  A  and  J^(o^  j(U))  are  linear  on 
every  edge  of  A  . 


Putting  i’  ~  n  ®J  (6.8)  gives 
P  P  P 


(7.11)  B  (n  ,v)  =  B  (e,v)  V  v£H  (<J>  ) 
P  P  P  P  P 


Using  Lemma  7.2  and  proceeding  in  the  standard  fashion  it  is  easily  seen 
that  n  is  the  unique  H  (<t>  )  function  satisfying  (7.11). 

p  p  p  }  ° 


LEMMA  7.3.  There  exists  a  constant  C  dependent  only  on  K  and  y  ,  A 
such  that 


(7.12)  1 | n 


P  1 
P  Hit  ) 


2 

1  C  f  l  ( | |g1+L1(U>) | | 
Li=l  H 


2 

+  l 

j=l 


Ar.)’] 


Proof .  The  set  of  possible  S^'s  is  finite.  Therefore  we  have  by  Sobolev 
imbedding  theorem 


(7.13)  | | v I 


I  o  l  c  I  i v  1 1  l 
H  (r.)  h1(*  ) 

P 


with  C  depending  only  on  K  .  (7.13)  with  the  Lemma  7.2  and  Schwarz's  Inequality 

leads  to  the  desired  result. 


LEMMA  7.4.  Assume  that  €  H  ( Z ,  p ,  c ,  ft)  then 


(7.14) 


|  >c£( i|g+L(U)||  +  Z  llixtf  (U))||  ) 

H  ■>  )  i=l  H°(4>  )  j  =  l  1  1J  H  (Ti) 


where  C  depends  only  on  K  ,  the  family  U  and  A,  y 


Proof . 

Suppose  that  the 

constant  C  does  not  exist. 

each  n 

=  1,2,3,  ... 

i) 

a  K-mesh 

ii) 

P[n]€R(P[n]) 

iii) 

i("ie"1<4p(„!> 

with  U 

^  being  bilinear  on 

iv) 

*lInl£L2(*p[n]> 

with 

II-n]!  !  1  c| |i[n]- 

1  H  (A)  1 

any  A  €5  with  C 

pin] 

average  of  on  A 

independent  of  n  and 

C  Pin] 

I  ft 

H°(A) 


such  that  the  unique  solution  n  €H  r  ,($  r  ,)  of 

pin]  Pin] 

2 

r  i  *-  .  r  i  r  1 


(7.15)  B 


+  j  <-'i<0i,j<U  »'Vj»-.,  In]  '  ^Vnl'W  • 

J-l  L2'  i 


satisfies 


(7.16)  |h[n]||  <  l  l  (||g|nl+  L  (Urn])|| 


H<«  ,  ,)  "  1-1 

P1"1 


“  (*  tn]’ 

P 


+  l  IU,(»1  >'",»ll>  0 

j=i  1  ,J  H°(rj  J) 

l 


Putting  r[nl  =  i[n]  +  L^U1"1) 
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£.!”!  =  J.(o.  .  (Uln!)) 
1,1  1  l.J 


then  without  loss  of  generality  we  can  assume 


(7.17)  Y 
i  =  l 


r . 

1  l 


11  o  -  1  ikin{M  0  - 

H  <%(„]>  H  "Y 


and 


(7.18)  |  r  [ n  1  |  .  <  C  |  |  r  ^  -  Mnl|  | 

1  V(A)  "  1  1  H°(A) 


VA  €S 


n] 


The  inequality  (7.18)  follows  from  the  assumption  that  g^ €  It  (Z ,p  ,e ,Q)  and  that 
L^(U^n^)  is  constant  on  A  .  From  (7.17)  we  have 


(7.19)  i|r' 


A-“[n]  |  | 

1  H°(A) 


1  \\*[ 


H°  (A) 


1  I  I  i 


<  1 


H  <*p> 


Using  Lemma  4.3  we  may  assume  that  <1> 


p(n) 


*  ’  S  [n]  =  S>  H  r„,(*  f^)  =  H 


Jn]  JnV 

P  P 


Further  by  Rellich's  lemma, (7. 18)  and  (7.19),  we  may  also  assume  that 


r .  ♦  in  H°(  i)  . 

t  i 


Now 


i  ,  j  i  s 


i.J 


yie Ids 


is  linear  on  each  edge  of  A €S  and  so  we  may  assume  also  that 
H°(l’)  where  T  is  the  union  of  all  edges  of  A€S.  (7.17) 


(7.20)  }'  (  !  |  r 


i  =  l 


I  o  +  l  N^’i  ill  o  )=  1 

H  (*)  .1-1  ’’  H  (f) 


(7.15),  (7.16)  and  (7.20)  leads  now  to  the  contradiction. 
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LEMMA  7.5.  Let  p€R(P)  be  such  that  =  (v^  >vp)  €  H (a)p)  (see (6. 7)) 


and  u  CQ.  for  some  k€Z  .  Then 
_ p  _ o _ 


I  h 


PUul 


2 
C  l 


H  (<S>  )  1=1 

P 


+  l  I  I 1(u>  I  I 


I  ®i  ‘  o  L  1  1  ~  ivi  \ ''  '  '  1  o 

h  (*  )  j=i  1  h  (r.) 


A ; 


where  restricted  to  any  AC$  is  the  average  value  of  g^  on  A  . 


Proof.  The  proof  will  be  by  contradiction.  Suppose  the  lemma  does  not 
hold.  Then  for  each  n  =  1,2,3,...  we  can  find  pfn^  ,  and 

gjn^  as  in  i) ,  iii) ,  and  iv)  of  the  proof  of  lemma  7.4  when  ii)  is  replaced 
by  ii') 


ii')  P^  €R(P  n  )  ,  (v  ^,v  jnj)€H(w  [n])»w^n]CQe  for  some  k€z 


and  (because  of  (7.10)  and  Bp(e,v)  =  0) 


(7.21) 


-<itnl  +  v  .j-J  >  o 

1  pln]  pLnJ  H°(4>  [n]) 

Pl 


,[n],N  _T-1 


=  0 


pc  1  pl  J  h  (r  ) 

i  Jn 1 
i » P 


.[n] 


and  with  as  in  (7.15)  and 


(7.22)  |  I  n 


[n] 


>  n 


£[n]  Altn) 


H  (%[n]>  Ll=1 


I  (||ginJ  -  8 


H  (*_ r_1) 


P  [  n  1 


+  i  iiv°i  j<^tn,))  n  o 

j=i  iJ  "hV  )J 


i.P 


[n]' 
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As  in  the  proof  of  the  previous  theorem,  by  use  of  lemma  4 . 3  we  can  assume 
that 


a)  S  r  i=  S  ,  <J>  r  ,  =  4>  ,  T.  ,  ,=  P.  ,  H  r  ,  =  H  are  independent 

fn]  pfn]  i.pfn]  l  p|n] 

P 


of  n 


f$) 


(1,  r  C  Q  for  some  k£z  where  k  is  independent  of  n  (because 


•pfnl  <9 
ZQ  is  finite) 


Now  a),  6)  yield  that  also  Bp[nj('»')  =  B(*,*)  does  not  depend  on  n 


')  h 


fn) 


i  1 

H  (4>) 


=  1 


Let  us  set  now 


d"1-  »'nl  +  Lja.'"1) 


,  .  j.(o.  .(iln)>) 

1,J  J-  1,1 


Now  by  (7.21)  and  y  we  have 


,[n]  _  A--  [  n  i 


1  H°(A) 


-*■  0  ,  VA€  4> 


(7.23) 


iu.  .it  -> o 

X’J  h  (r.) 


t  l 

f 


\ 

l  — 


It  is  easily  demonstrated  that 


(7.24)  l  | |L  (U[n])  -  %  (U(n])  M  1  l  IU  || 

i=l  H°(4>)  i  ,  j*=l  H  (T.) 


and  so  we  must  have 


(7.25)  j  |l_.  (l'|ni)  -  VL.  (l'lnl)  ||  ,0 

V(t) 


By  (7.22),  using  lemma  7.1  Rellich's  lemma  and  the  fact  that  S  is  finite, 
we  mav  assume  that 


(7.26)  -»  r\  in  H°(t) 


(7.27)  ijn|  -  g.  in  1!°(;) 


Because  of  lemma  7.1  wo  see  that  g. 


also  in  H^^(i)  where 


H  (P)  is  the  usual  space  of  fractional  derivatives.  Now  by  (7.23)  again  we 
see  that  r^  and  g;  must  be  constant  on  each  A€S  .  Further  by  lemma  4.3(i) 
and  the  fact  that  g.€H^*(4>)  we  see  that  g^.  is  in  fact  constant  over  all  of 


By  (7.26)  and  (7.27),  it  follows  that  1,^(1’^)  converges  in  H°(|l>)  and  by 
(7.25)  it  must  converge  to  a  constant  on  ■>  .  This  follows  that  r^  is  a  constant 


By  (7.21)  in  the  limit  we  see  that 


/v  .-1 

(r.,v.  °J  .  )  =0 

1  Pf°l  P["l  H"f*) 


Bv  lemma  1.5  we  have  v  .  . >  0  and  v  r  ,>  0  on  a  set  of  positive  measure. 

Pin]-  Pin] 

Or 

Thus  we  conclude  that  r^  =  0  . 
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Returning  to  7.15,  we  see  that 


In} 


in  H1  ( 4  ) 


which  contradicts  Y  • 


Now  we  return  to  the  question  of  riie  u-posteriori  estimate.  Oiven  a 
mesh  V  and  the  corresponding  finite  element  solution  l'  ns  in  (6.5), 
we  will  associate  to  every  t-€0  the  error  indicator  n(A)  defined  by 


2  2 

(7.28)  n2(A)  =  1  {diam  A  £  |  |  J.cj  .  .  (U)  |  |  “  +  (diam  A) 2 | | g ,+L . (U) | | 2  } 

i=l  j  =  l  1  1J  ll°(a.A)  1  1  H°  (A) 


where 


? 

a-u 

Vu)  " 


1  ~2 

au. 


L2(U)  =  (A+u) 


au  au  au.  au. 

o.  .  =  x(~-  +  -_— )  +  p[-— -  +  -~J ] 

l,j  ax,  3x_  Ox.  3x. 

12  ji 


and  (J2)  indicates  the  jump  across  the  vertical  edges,  3^A  (horizontal 

edges,  32A  )  of  A  .  We  use  the  same  convention  as  before,  i.e.  if  the  edge 
rein  is  such  that  v .  £  H .  -*■  v .  =  0  on  T  then  the  lump  J.  is  not  taken 
into  consideration.  If  it  is  not  the  case,  we  take  the  jump  J.(d.  .)  =  a.  . 

1  1,1 

The  estimator  E  then  is  defined 


E2  =  l  r,“(A) 

A  €0 

THEOREM  7.6.  The  estimator  E  is  an  upper  and  lower  estimator,  i.e.  there 
exists  c^,  c?  dependent  only  on  K,  X,  11,  U  and  f  but  not  u  and  V  ,  such 
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that  the  error  of  the  finite  element  solution  of  the  problem  P(HQ,w,g) 
w€[M(P)]^  ,  g  (Ztp  ,t  ,fl)  satisfies  the  estimate 

(7.30)  c  1  <  | )e| |  <  c,E 

1  1 _ 


Before  proving  the  theorem  let  us  prove  two  simple  lemma. 

LEMMA  7.7.  Suppose  F  is  an  edge  of  some  A£P  .  Then 

a)  Either  (i)  rC3ft  , 

or  (ii)  3  Pf  R(P) 

such  that  ACo)  but  F<I3u) 

P  ~  P 

g)  F  fl  A  ’  ^  0  for  at  most  x  elements  A'6D  where  x  depends  only 
on  K  . 


Proof . 

a)  By  lemma  3.5  we  have  Ev^  =  1  ,  so  certainly  for  some  P  €R(P) 

v  >  0  at  the  midpoint  of  T  .  Since  v  is  a  non  negative  and  linear  on  T, 

P  P 

then  Vp  >  0  on  f  except  perhaps  at  an  endpoint.  The  result  is  then  immediate 

upon  noting  that  v  is  continuous  and  inducing  lemma  3.6. 

P 

0)  The  number  of  such  elements  is  clearly  bounded  by  the  number  with  non 
empty  intersection  with  A  .  Lemma  4.1  iii)  then  gives  the  result 

LEMMA  7.8.  Suppose  a^,...,am  are  non  negative  real  numbers  and  let  a 
be  a  mapping  from  {l,...,n}  -»  {l,...,m)  such  that 
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i)  a  is  onto 

ii)  ]  H  '  0  such  that  Vl  £  j  1  m  the  set  {i|l  <_i  j£_n,  a(i)  =  j  }  has 
at  most  M  elements  then 


m 

l 

j  =  l 


a. 

J 


n 

Ia«.  ( i ) 


m 

1  M  l 
j  =  l 


a  . 
3 


The  lemma  is  obvious. 


Proof  of  the  Theorem  7.6.  Using  a  scaling  argument  applied  to  lemma  7.3  and 
7.4  we  obtain 


(7.31) 


» <v 


L  1  [ (diam) 2 | |g  +L  (U) 
A€D  i  =  l  11 

AC-o 

P 


H°(A) 


+  diam  A  l  ||J?(o  (U))||2  ] 

.1  =  1  1  1J  H°(d.A) 


where  J°  indicates  that  we  use  the  jump  across  edges 
but  within  fl  equal  to  zero.  Adding  over  all  P€'R(P) 
and  lemma  4.1,  we  are  in  a  situation  covered  by  lemma  7 
Before  formulating  the  next  theorem,  we  prove 


A  which  are  on  3w 


P 

and  using  lemma  7.8 
7  and  the  result  follows. 


LEMMA  7.10.  Let'  R*(V) 


{P  €R(P)  j  v 

P 


0  on  U  3Q^  }  .  Then 

k€-  7  9 

o 


(7.32)  U  *i.  »  Q 

P€R*(V)  p 
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Proof.  The  lemma  follows  from  lemma  3.2  applied  to  each  separately 

■ -  d 

which  shows  that  any  A€V  has  at  least  one  vertex  P  say,  being  a  regular  node 

k  k 

located  in  the  interior  of  a  .  Clearly  the  corresponding  ,  and  so 

P€R*(V)  . 

Let  us  introduce  now  another  error  indicator  n(A)  , 


(7.33)  n2(A)  =  J  {diam  A  ^  | |J  o  (U) | | 2  +  (diam  A) 2 1  | g.-AgJ  1 2  > 

n  i  .  •»  *  9  J  _  »  1  1  ..O  /  .  v 


«.=  !  j=l  1  1,J  H°OiA) 


H  (A) 


where  g^  is  the  average  value  of  g^  on  and  the  corresponding  error 


estimator 


(7.34)  E  =  l  n  (A) 
A€V  1 


THEOREM  7.11.  The  estimator  (7.34)  is  also  an  upper  and  lower  estimator. 


Proof.  Let  p€R*(V)  .  Then  using  lemma  7.5  and  7.4,  we  see  that  for  any 


(7.35) 


I  (diam  A)2| |g  +L  (U) | |2  <  f  ( (diam  A) 2 j  | g  -Ag  | | 

i-1  H°(A)  A6w  i=l  1  1  H°(A) 

P 


+  diam  A  f  ] |J°o  (U) |  | 2  ) 

j  =  l  1  1,J  H°O.A) 


We  have  shown  in  lemma  4.1  that  any  A  is  contained  in  not  more  than  M(K) 
stars  ii'p  .  Lemma  7.10  shows  that  every  A  is  contained  at  least  in  one  o> 
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P  fc  R* (P)  and  therefore  summing  (7.33)  over  all  A  CD  we  see  immediately  that 

-  ■v  m  —  a. 

K  <  CK  which  yields  that  E  is  an  upper  estimator.  Because  E  <  E  E  is  clearly 

a  lower  estimator. 


Lemma  7.12.  There  exist  CL  >  0  ,  i  =  1,2  such  that 


(7.36)  C  ||J  n  .(1011  o 

H  '(). A) 


(A)m2 


-t 

<  diam  A  )  fJ.o  •  (U)  (x^'' )  ]  < 

Jt=l 


c2hVi,j<l'dl'. 


H  (?.A) 

l 


(/\) 

where  x?  are  tlie  vertices  of  A 


J.  r, 

1 


Proof. 


.  .(!.')  is 

i ,  J 


The  inequality  (7.36)  follows  immediately  from  the  fact  that 
linear  on  every  edge  of  .A  . 


Lemma  7.12  allows  us  to  introduce  another  error  indicator, 

2  2  4 

(7.37)  n2(7)  =  )  {(diam  A)“  )  j  [J.o  ,(l!)(x^)]"  +  (diam  A) 2 1  |  g . -Ag .  1  |  2  } 

•■■■  1  L  j  =  l  f.=  l  1  1,1  *  1  1  H°(A)J 

and  error  estimator 

(7.38)  E2  =  }'  n-'(A) 

At"  ' 

and  wo  have  from  lemma  7.12  and  Theorem  7.11, 


THEOREM  7.13.  The  estimator  (7.38)  is  an  upper  and  lower  estimator. 
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8.  THE  ASYMPTOTIC  ESTIMATE 


We  have  shown  in  the  previous  section  that  the  estimator  E  introduced  in 
(7.35)  and  (7.34)  is  under  the  proper  assumptions  simultaneously  an  upper  and 
lower  estimator. 

In  this  section  we  will  analyze  an 'estimator  E  ,  which  will  be  equivalent 
with  E  i.e.  there  exist  and  so  that 


(8.1)  C2E  £  E  <  CLE  , 


and  will  be  asymptotically  exact  for  the  energy  norm  |||’|||  introduced  in  (6.3). 
We  shall  say  that  an  estimator  E  is  asymptotically  exact  with  respect  to  the 
energy  norm  if 


(8.2) 


i  •+  1  as 


To  show  (8.2)  we  have  to  make  various  assumptions  about  the  solution  u  and 
the  meshes  in  addition  to  the  assumption  that  the  mesh  is  a  K-mesh. 

Suppose  V(Q)  is  a  mesh  and  let  P'CP  satisfy 

i)  If  ACC'  ,  diam  A  =  h  ,  h  >  0 

ii)  if  A  CD'  then  all  vertices  of  A  are  proper  nodes. 

If  (i),  and  (ii)  hold  then  we  shall  say  that  V'  is  uniform.  If  only 
(1)  holds  then  we  will  say  that  V ’  has  uniform  size. 

Suppose  now  that  A  C  B  C  n  then  we  shall  write  A  <  B  if 
P(A,B)  =  dist(A,lR2-B)  >  0  . 

Let  PQC D(n)  then  we  will  write 
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(  l'  )  r  I  n  t  1 

■€D 


.lnd  (( )  )  M  iiiiln  i  ■  >;  el  i-ments  rout  a  ined  in  V  .  We  will  write  V"  <  V '  if 

o  a 


(V")  . (P‘ ) 


Finally  denote  bv  M(P;’  )  the  set  of  functions  of  (M(P)]  restricted 

O 

to  and  M(  P,  )  C  M(p,  ’  )  the  set  of  functions  which  are  zero  at  . 

We  shall  Make  use  of  the  following  version  of  Theorem  5.2  of  111]  . 


THEOREM  8.J_.  Let  0  ’CP(i;)  he  a  uniform  mesh.  Assume  that  the  bilinear 

form  B  defined  in _ (6.1)  has  constant  coefficients  on  n(P')  and  let  V” <  V  . 

Then  if  \  =  (1  ^ ,  r-,)  €  I  H^  (<  (P' )  ]  "  satisfies 


(8.3)  i b '  -  A  >  j:’ 

i=l  H 1  (:..(£?'>  > 


for  a  L 1  l^€f.l(V,u(D’))  then 


U-~  x| 


( 8  •  A )  I  f ■»  !  ,  ll  C  A  f  inf  i  L  =  '  •<  !  I  +  n(fj  ’  )  )  I  ^  -i  Xl  0 

'  H  (fUP")  'l  =  l\y€M(P  >}  H° (ft(P ' ) ) 


(:i(V'r)  ,’  .(0r) )  1  j 


II  (:  ;(P' ) ) 


with  C  dependent  only  on  the  b  i  I  inear  _f  orm  (6.1)  (and  not  on  fl(P")  ,  Q(V)  )  . 

Suppose  now  that  P'CP  has  uniform  mesh  size  and  let  P"  <  P'  .  Then 
we  can  defim  a  d  i  I  ferenco  operator  I.  wliich  maps  H^(U(P'))  into  H  (fl(P")) 
so  t  ll<)  t 


(T.u)(x>  =  2h  |  ii  (  x+ho  .  ) -u  (x-lie  •)  | 

i  l  1 
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with  h  being  the  diameter  of  the  elements  of  £)'  and  e^  the  coordinate 
unit  vector. 

We  have  then  the  following  theorem  which  follows  easily  from  Theorem 
6.2  of  [  ]  . 


THEOREM  8.2 .  S uppose  that  u^.  and  U.  j  =  1,2  are  the  components  of  the 
exact  respective  finite  element  solution  of  the  problem  P(HQ,w,g)  as  defined  in 
Section  6. 

Assume  further  that  V"  <  V,  P'Cp  are  as  above,  that  the  bilinear  form 

3 

B  of  (6.1)  has  constant  coefficients  on  Q(D')  and  Uj  €  H  (D(D’))  ,  j  =  1,2  . 

Then  if  V" '  <  V"  it  follows  that  for  D.  =  -r-^—  , 

l  dx.  we  have 

l 


(8.5) 


|DiurTiuj 


H°(S1(P"')) 


2 

C  I 

k=l 


(h 


+  p  (n(p,"),n(p")) 


(Q(P’)) 


) 


Denote  now  for  any  square  Q  C  IR  by  (Q)  the  space  of  all  functions  which 
are  of  the  form  ax^  on  Q  (when  referred  to  an  origin  placed  at  the  center 
of  A  )  . 

Now  we  prove 


THEOREM  8.3.  Using  the  same  notation  as  in  Theorem  8.2  then  for  each 

A€P"'  there  exists  <t>^  .  €H.(A)  i  j  =  12  such  that 
_ hi _ i _ _ _ : _ 


)  !  |D 


(u  -U  )  -  l  <p  ||  <  Cp_1(fi(P’"),fl(P"))  l 

3  3  3,i  H  (n(p,M»  k-1 


W\K\]  ■ 

K  H‘ 


(n(P"’ )) 


i  I  uk"\  1 1  o  } 


H  (fi(PM)) 


8.6 


6A 


Proof.  For  any  w  0!(P,  d(D") )  wo  have 


(8.7;  B(l:  -  u,w)  =  B(u-!:u,w)  =  )  B,(u-i>u,w) 

A€P"  A 


where  ilu€ 'I(P)  is  such  that  (  !u)  (x)  =  u(x)  for  any  x€R(V),  and  is  the 

restriction  of  B  to  A  . 

On  each  !: \€V "  ,  B.(u-.'iu,w)  is  the  sum  of  terms  of  the  form 


[  'W  . 

1  U  — — -  (u.-P.U.)  v~ 

I  -»x.,  1  .1  -X 


for  =  1,2  ;;nd  some  constant  ■<  .  Consider  first  the  terms 


I  on  . 

J  Vx,  (V  ’’ 

J  K 


iw  . 

)  —  1  (lx 
ax, 


r  f  is  ensv  fn  that  for  a  quadratic  function  the  above  expression  is 

equal,  to  zero.  Therefore  by  standard  argument  we  get 


>w . 


(8.8)  i  -  (u.-u.)  ’x  ’  Ch"|u.  Iw.  . 

1  ’  ^  1  V(A)  1  V(A) 


with  C  independent  of  and  u,w 

hot  us  consider  now  the  term 


aw  . 

1 


(8.9)  ;  —  ( ii  .  —  ’ii .  )  -----  dx,  k  ^ 

,  .ix  1  1  ax. 


Obviously  we  need  only  analyse  the  ease  k  =  1  ,  ?  =  2 


I'sing  the  notation  as  shown  in  figure  8.1  we  integrate  by  parts  in  (8.9) 


Figure  8.1.  The  notation  of  an  element. 


(8.10) 
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SW  . 

Ren  I  i a  i  m>,  tii.it  (u  — i;u  .  )  is  a  ivnt  inuous  f-.inot  ion  on  1  (D")  and  •- —  Is 

.1  J 

<w  .  1 

continuous  on  all  horizontal  oip.es  of  V"  and  =0  on  the  horizontal  edges 

'X  I 

of  ... 1 P")  wo  got  by  adding  (8.10)  over  all  €  V" 


(8.11) 


.  >') 


(u.-fiu.)  -  dx;  ■"  Ch"f<!.  ;  |w,J 

xk  1  3  xi  “  1  u\n(V"))  1  H  (cjTV")) 


for  k  i  i 


(8.8),  (8.1|)  and  (8.7)  viold 


•  8 .  |  J  )  lid'-  u.w1  <  it '  (  ii.  ,  )  (  /  '  w .  I  ) 

j-  i  '  m 'i  ft;” i i  1=1  1  ir(o(P")) 


Applying  now  t  hoor-u:  >8.1)  wo  got  i  rcni  (8.111) 


(8.1  i )  .I/.-’  <  i .  !  -l 

1  1  11  (  CD"’  )  j 


1  +  (o(D")  nCD'"  )  ^3 i_'7'u i  I  o 

H  (::(P"))  1  ’ '  KU  ’  J  3  H  (8(0" 


(  (O' 


•y  ’  {h" !  'uj  i 1  )  +  I  i!’Vu  J  i  o  j 

Y; - 1 L  1  II  '(:'(£>”))  1  1  H  (P(h’"))J 


whore  wo  have  used  t  I,.  Out  that 


(8.14)  ji  -’a.  1  1  el  \  V~M„  i]  +  j  |  u.-l!u  .  1  I 

1  1  !i<r(  ■iV"))  1  1  u°(i.(V"))  3  •’  H°(;?(P") ) 


-  n,2l|..|!  , 

h  (,'(P")>  ir(w(P”)) 


a  rtf 
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Let  us  consider  further  D.(uj-riu.)  on  A  £0"'  .  Refering  to  Figure  8.2 


Figure  8.2.  The  notation  of  an  element. 


2  A 

for  any  w£H  (A)  we  choose  4>^  =  J^w€H^(A)  so  that 


(Z(l))  =  w(Z(l)) 


where  Z^^  are  the  midpoints  of  the  particular  sides  as  shown  in  Figure  8.2. 
For  u.  =  g  a  quardatic  function  ,  we  can  easily  check  that 


And 


o^g-ng)  =  JiDig  . 


so  by  the  usual  arguments  we  have  that  for  some  ,€H. (A)  , 

3  * 1  1 


(8.15)  ||D  (u  -flu  )-♦*  Jl  <Ch2|u  j 

1  3  3  J>1  H  (A)  3  HJ(A) 


Summing  over  all  &€V" 


we  obtain 


(8.16) 


1).  (u  .-[Til  .  ) 
1  J  1 


i ,  i 


it 

u  ('. 


Clf 


{p"' )) 


;> 


Combining  now  (3,1  6 )  and  (8.13)  wo  t  the  desired  result. 


KKMMA^  6 .  A .  L'sins  tin-  sar.it  notation  ajs  in  Theorem  8.3  we  get  for  k 


Coin 


i  a  .  -  l  .  I  t  | 

1  1  !i°(  av" )) 


h  |  j  u  .  i 


H3(ii(D")) 


( V" )  i 


h" 


'n3, 


CD”)) 


|urujii 


W!1  ‘TC 


f  ( V"' 


( 


•V* 


>  ) 


P rnoi  ,  W r i  ic 


(■  (O'"  ) 


(p"> )  < 


j= i 


1!  \  (V"  )  ) 


I 


ir  (a  (£>")) 


Now  we  have  ior  anv  \  €  D" 


(8.17)  fn.i  a. ->■.),[)  .(d.-c.))  =  (;•.  ,+n,(u.-r.)  - 

i  I  K  1  )  .  .  ]  .  ■  1  1 


II  (\) 


I  .  ■  v  -I  1 


+  °k(,V 


,  ,  >  -  ( ,  )  +  (n.(u.-u.)  -  /  „ , <&’!'  ,  ) 

’’h  n°c,)  ’ k  h"(a)  !  1  »  J>?  1  ,k  ll°(A) 


4  (u.-v.)-t:  .  )  „  +  „d  „(u, -(-)-*'  ,) 


±A 


1 


69 


It  is  easy  to  see  that  for  P  ^  k  we  have  (<{k  ,)  =  0  .  Summing  (8.17) 

J*4  i,h  Ho(A) 

over  all  a€V"  and  using  Theorem  8.3  we  get 


(8.18)  |  (D  (u  -U  ),D  (u  -U  ))| 

'  J  J  k  1  1  H°(Q(V"')) 


<c«»2-n.(  l  l  |*‘  I2  )1/2) 

i,j  =  l  A  €.  V'"  1,3  H°  (A) 


Using  once  more  Theorem  8.3,  we  get 


(8.19) 


v  A 


J.1'  U° 


tT(fi(P"'  )) 


<  C[<H-|  |D  (u  -U  )  I  I 


1  J  J  u° 


Hu(n(P'” )) 


Using  now  Theorem  8.1  for  £  =  u-U  we  have 


(8.20)  ||D  (u-U  )|] 

1  J  3  H°(n(p'”)) 


c. 

£  c(  l  (h| |u  I  I  + 

j=i  3  ir«i(P")) 


3  3  H°(fi(P")) 


and  combining  now  (8.18),  (8.19),  and  (8.20),  we  obtain  the  desired  result. 


LEMMA  8.5.  Let  S  be  the  square  as  in  Figure  8.1  with  side  length  h 
Then  for  any 


f  =  ax^  +  bx£  +  cx.^2  +  8 


we  have 


(8.21)  |  |  ax.  |  j  2  -  h2(f2(N+)+f2(lO+f2(S~)+f2(S+))| 

1  H°(S)  W 


J 

<  M  ||bx9+cx.x9+d  |  I 

H°(S) 


with  M  independent  of  h,  a,  b,  c,  d. 


71 


Proof.  We  have 


(8.24)  D.(u.-U.)  =  (D.u.-ThU.)  +  (ThU.-D.U.) 

i  J  J  ijij  i.]ij 


By  Theorem  8.2  we  get 


(8.25)  1  |D  u  -  tJu  ||  IC[  l  (h2||u  ||  3  + 

J  1J  H°(U(P"'  ) )  k=l  HJ(fi(P')) 


+  [p(fl(P"’ ),i7.(0"))f1|  |u  -U  |  |  n  ] 

R  R  H°(fi(P')) 


Theorem  8.3  yields 


(8.26)  ||d(u-U)-  l  /  .|| 

J  j  /\€P"«  h  (n(P"' )) 


<  do"1(fl(0,n),n(P"»  l 


k=l 


{h2||\ll  3  +  I  I  Vuk!  I  o  >1  =  E 

hj(\.!(P")  k  H°(n(P"))  1 


Combining  (8.24),  (8.25)  and  (8.26)  we  get 


(8.27) 


j  |t!>  -D.U.-  >  *711 

1  J  )  1  a dv”  i  »x  u° 


<  C  E, 


•€P"'  1,1  Hw(fi(0"' )) 


Simple  computation  shows  that  T.U^-D^  is  bilinear  on  each  A€0"'  .  For  sake 

of  definitions  and  without  loss  of  any  generality  let  us  suppose  that  i  =  1. 

Consider  now  T^Uj-D^Uj  .  Using  notation  shown  in  Figure  8.3  we  have 
for  F,  =  (fn.^2)€rE 


r 
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N'  NT 


Figure  8.3.  The  scheme  of  the  notation. 


(8.28) 


(tijV-D.lhHO  =  2h  1tUj(f.1+h,C2)  -  IKC^-h,^)]  -  h_1tU:j(0-U:.(F1-h,C2)] 


1 

2 


For  £6ru  we  get  analogously 

(8.29)  (tV-D  U  )(C)  =  -  \  J.  AO 

h 

Let  now  <J>  be  the  L„  projection  of  (T  U.-D  U.)  onto  H  (A)  .  Then  from 

^  *  J  J  -I* 

(8.27)  we  get 


(8.30) 


lTiYDiV 


h  (n(P"' ) ) 


:  I 

&€V"’ 


iK^2o  +o(ei> 

H°(A) 


Now  it  is  easy  to  see  that  for  any  bilinear  function  f  on  A  ,  f  *  ax^  + 
bx2  +  cXjXj  +  d  (with  coordinate  origin  in  the  center  of  A  )  the  L2  projec¬ 
tion  of  f  onto  H^(A)  is  exactly  ax^  and  so  by  lemma  8.5,  (8.28)  and  (8.29) 
we  obtain 


(8.31)  l  ||$  ||  =  fg  l  P  (N_)  +  J  (N+)  +  J  (S  )  +  J 

/\£P"<  1 > J  u°/*\  48  a^D'"  *>3  1 


H  (A) 


A€Z? 


0( I jT^U  -D  u  -  l  ||  )  . 

J  33  A€P"'  1,3  H  (n(£>"'  )) 


using  now  (8.27)  and  (8.30)  we  get 


(8.32)  ||tV-DU  || 

13  3  H0(n(P"' )) 


<l>?  .(0"’ )  +  0(F.J) 
A  *  J  A 


Hence  we  have 


(8.35)  | (D  (u  -U  ) | | 

1  3  3  H°  (;’«?' 


I  I  (T!?U.-D1U.)  +  (D1u.-T111U.)  ]  I 
>)  1  3  1  3  1  3  1  J  H°  (8(0'"  ) 


I  T^U .-D,U .  '  |2 


13  1  3h<W)) 


+  Mv  -Tiu  ll2o  + 

13  L  3  H  (.HO"1  ) 


+  «| |D  u  -T  U  | |  I|T"U  -D  U  || 

13  33  H°(ij(P"' ) )  3  3  H0(n(P"' )) 


where  -2  <  n  <  2 

Further  using  (8.25)  and  (8.20) 


(8.34)  ||tV-DU||  <||d(u-U)M  + 

1  3  33  H°(n(0"))  3  3  3  h  (n(P'" )) 


+  T 


lUl“Dlull  I  o  -  C(E1  +  l  hMUi  I  I  2 

13  3  H°(fi(P"))  3  j=l  3  n 


(n(V)) 


+  P~1(n(V" 


|  u . -H . I | 

3  3  H°(S?tP')) 


,(S+))  + 

9  J 


>,n07")) 


Using  now  (8.32),  (8.33),  (8.34)  and  (8.25),  we  have 
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(8.35) 


IVWll'o  =  i^”')  +  0(eJ  +  E  (  l  h||u  ||  + 

1  J  3  H°(P")  1>J  1  1  j=l  J  H3(J2(P,)> 


+  p"1(fi(P",),n(r>)||u.-u  ||  ) 


3  3  H°(fl(P’)) 


and  (8.35)  yields  almost  immediately  the  theorem. 

In  section  6  and  7  we  introduced  various  error  indicators.  See  e.g. 

theorems  7.11  and  7.12.  It  is  obvious  that  we  can  have  many  equivalent  error 

indicators  which  would  be  simultaneously  upper  and  lower  ones.  Theorem  8.6 

enables  us  to  design  a  special  one,  optimally  suited  to  our  purpose. 

In  (6.1)  we  introduced  the  basic  bilinear  form.  Obviously  we  can  write 
T 

B(u,v)  =  [Dv]A[Du]  where 


[Du]  -  [D1u1,D1u2,D2u1,D2u2] 


and  analogously  Dv 

The  matrix  A  has  then  the  form. 


- 

- 

— 

\+2u 

0 

0 

A 

0 

u 

u 

0 

A11 

A12 

0 

u 

u 

0 

A 

A+2p 

A21 

CM 

CM 

< 

_ 

_ 

_ 

J 

Assume  now  that  A€P"'  as  in  theorem  8.2.  Then  define 


Q  =  j^1{fJl,l’Jl,2lAll[-7l,l’Jl,2^  [aj]  + 


[aJ]} 
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where  a.  j  =  1,2, 3, A 
Further  let 


are  the  four  vertices  of  the  element 


A  . 


f^(A) 


i 


i=l 


2 


H°  (A) 


and  then  define  the  error  indicator  (which  obviously  is  simultaneously  an  upper 
and  lower  one) 


(8.36)  n2(A)  =  |g  a2 (A)  +  h2-Y82(A) 


where  y>0  is  a  constant  which  will  be  determined  later.  We  mention  here  only 
that  for  smooth  g  we  have  J  3 ( A)  |  £  Ch  . 

Now  we  have 


THEOREM  8.7.  Let  the  assumptions  of  Theorem  8.6  be  satisfied.  Then 


(8.37)  HI u-U | 


.1(17"  ) 


=  I  a2(A)  +  C[  l  { (p  1h3+p"2h4)f|u1  ,t2 


Ac'?" 


2 

l 

k=l 


k  1  1  ^ 

HJ(fi(0')) 


+  P~2  I  I  e  |  |  2  +  (p  1h+p  2h2)  |  |  u,  |  |  .  I  I  e  |  | 

H°  *  HJ(Q(P’))  H°'«' 


+  {  l  hl |u. | |  3 

j=l  J  HJ(n(P')) 


+  p  1 1  I  e  |  |  1  <  1  h'  ,  |U  |  |  + 

H°(n(P'))  j=l  J  HJ(fi(P')) 


+  e 


} 


Hu(ft(P’)) 


Proof.  We  have  for  u  -  U  =  e  = 


l  Mill 


2 

«(P'M ) 
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j  \  (A+2ti)  ((Die])2+(D9e2)2)  +  u[  (U2ls1)2+(D]e2)2]  +  ZplD^D^]  + 

fid?’") 

+  2\ [D^e^D^e^ ) }dx 


Using  Lemma  8.4  and  Theorem  8.6  we  get  the  desired  result. 


So  far  we  have  defined  the  error  indicator  n(A)  by  (8.36)  for  all  A€  V”  . 


We  will  extend  this  definition  to  all  elements  A  with  A€Q  (see  sections 

u 

2  and  6)  where  the  coefficients  are  constant  and  A  D 3Q  =  0  .  If  A  has  an 

U 

edge  on  3Q  then  we  shall  define  the  indicator  in  terms  of  the  jumps  in  o.  .  , 
u  i » j 

so  that  we  obtain  an  upper  and  lower  estimator.  The  detailed  formulation  arid 
extensive  numerical  experience  will  be  addressed  in  a  forthcoming  paper. 

Here  we  will  assume  only 

a)  The  indicator  has  the  form  (8.36)  for  any  A  € ,  Afl3Q^  =  0  (y  >  0) 
fixed  but  arbitrary. 

b)  All  the  meshes  on  K-meshes  and  the  indicator  leads  to  a  simultaneous 
upper  and  lower  estimator  (under  the  assumptions  spelled  out  in  Section  7  ). 

We  will  place  some  additional  assumptions  on  the  solution  u  and  the  meshes 
used,  as  listed  below.  Then  we  shall  prove  that  the  error  estimator  is  asymp¬ 
totically  correct.  By  this  we  mean  that 


as  ||  |  e  |  |!  -*■  0  . 

We  will  assume  that  there  is  a  sequence  of  K-meshes  h(P^) 

h  -*■  0  with  K  independent  of  h  . 


h. 


:)  The  mesh  is  equilibrated  in  the  sense  that  max  n(A)/  min  n(A)  C 


with  C  independent  of  h  . 

d)  !|)e|!|  ;•  Ch  with  C  >  0  independent  of  h  . 

e)  There  exists  s  >  0  and  C  ,  both  independent  of  h  ,  such  that 


MIL  iCh^Mel 

H°(Q) 


f)  For  each  e  >  0  ,  li  >  0  ,  there  exists  meshes  V'  .  >  V"  ,  >  V 

h ,e ,  l  h , e ,i  h , e , 

i  =  l,2...m(h,e)  such  that 

i)  the  bilinear  form  has  constant  coefficients  on  ST ( P ’  .) 

h,c  ,i 

ii)  V"  .  are  uniform  and  V'  .  are  of  uniform  size 

h,£  ,i  h  ,£  ,i 

iii)  n(£>"'  .  )fl  .)  =  0  for  all  i  4  j  . 

n  » :  - ,  l  n ,  c  ,  1 

iv)  p(Q(P"  .),  .))  >  Ch6  with  C  and  9  independent  of 

h ,  t  ,  l  h ,  e  ,  i  — 

h  and  £  ,  0  <  9  <  S 

m(h,t) 

g)  Denote  Rl  =  U  S2(P'  . )  and  analogously  R"  and  RT’  and 

n ,  e  i  —  ]  n  ,£  ,1  u ,  c  rijc. 

we  shall  assume  that 

0  I  M  |  ,  1  X(e) 

h  ,£- 

ii)  |  | g |  |  ,  i  X (e) 

>•  (Rh,c> 


m(f  ,h)  6(h) 

h)  Let  S’"  =  9.  -  U  SKPL"  •)  and  S"'  be  the  6  neighborhoou 

h,£  h,e,i  h ,  t 

e=  l 

of  S"'  with  *  =  h°  0  <  9  S  . 

h ,  £ 

Denote  now 


&  .  =  ( A  I A  fl  S"’  ,Uh)  t  0}  , 
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ancf  we  will  assume  that 


o(i)  )  <-  fo(P  ) 

h,t  —  h 


where  o(P^)  respective  a(ft^  f)  denotes  the  number  of  elements  contained  in  P^ 


respective  ft, 

h ,  f. 

Now  we  have 


THEOREM  8.8.  Suppose  that  the  assumptions  listed  above  are  satisfied,  then 


e  >0  such  that  for  e  <  e  and  h  <  h  (c) 
o _ _ o _ —  o 


<  Mr 


where  M  is  independent  of  h  and  t: 


Proof.  For  every  Pj^'  ^  ^  we  can  write  by  theorem  8.7  and  property  f)  and 


g) 


m(h,r)  0 

?  ^e^(P'"  .) 

i=l  h,e,i 


l  n  (A)  +  z  -  2 


'  €  U  V"h  F  i 

i  =  l 


where 


3-O..4-20,„2,.,  .  L-20+2s|i|e(i|2+  ^1+8-9^24,-29,  x(e)  ,|j .  |||  q  + 


|z|<  c[(hJ-v+hH  "u)x"(r)  +  h 


1  (h1+S-6X(r)  +  hZS-Z0|||e|||o)(x(F)h"-S+|||e|||)] 


2s-26 


2-s , 
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and 


|  z|  = 


■yh 


m 

A€U  V" 
i=l 


J  I !gi'gi 

i  =  l 


h  ,  r  ,  i 


2 

H°  ( A) 


Taking  into  consideration  the  property  d)  we  get 


z  <  C 


(h1-fVh2-2s)X2(e)  +h2(s-9)  + 


(hs-e+hi+s-20 


/ul-e  2,  ,  ^  ,  l+s-20  ,  N  ,  us-6  ,  ,  .  ,  2s-20, 
(h  X  (c)  +  h  X  (  e  )  +  h  x(e)  +  h  ] 


We  have 


2 

H°(A) 


< 


Ch2|  |g. |  i 


2 

HL(A) 


and  therefore 


ch  1  I  I  ) 

i=l  h  ,c  ,  i 


Ch 


x(e) 


Hence 


\z\  +  \z\± 


,2 

'S2 


for  h  <  h,.(f) 


)  x  ( e )  + 


We  have  further 
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J.  n2(A) 
h 


and  by  property 


2 

n 

max 

2 

n  4 

min 


1  C cl  n2(&) 
A€D 

h 


when  using  property  c)  . 

So  we  have  for  sufficiently  small  e 


l  n2(A)  <  a(b  )n2  (A)  <  r.o(P  )n2  <  e  £  n2(A) 

£  h,t  max  -  h  max  -  in 

\€V.  Ath 


m 

U  n'" . 

r.  ,1/  h ,  c  ,  l 

i=l 


^  ■,  (A)  +  ,  l  n  (A) 


m 

,\€  U  V' 

i  =  1 


h  ,  r  ,  i 


MED. 

n 

.  m 

A f  U  O'"  . 

r  h,t,i 


h)  we  get 


(1-eM)  )'  n2(A) 


n2(A) 


m 

A€  U  P'" 
i=l 


h ,  c ,  i 


£  J  n2(A) 

”  A€Pu 
h 


Using  further  the  .act  that  the  estimator  is  an  upper  and  lower  estimator  we  get 


mOi  ,c)  2 

f » -  39 >  y  |  |  I e I  I  v  =  E‘4(l+B(h,e)r.) 

f  =  1  V  h  ,  f  ,  i  ’ 


(h,  )  j  ■  M  ,  independent  of  h  and  e 

j  I  v  we  have 


m(h 

) 

i  - 


I2 

^(P’" 


h ,  r. ,  i 


|2 

'S'" 


h  ,e 


— — 


i  2 

Let  us  analyze  now  |||e|||g. 


It  is  easv  to  see  that  there  exists  a  function 


such  that 


i )  0  <  ik  <  1  on 

h  — 

“>  iD\i  ;-s-fhTi  Ch‘° 


iii)  *h  =  1  on  S- 


ip,  =  0  on  fi  -  S"'  6(h) 
h  h,c 


Obviously  now 


(8.40)  I II  e  || |  s<„  £  B(eiJ)h,eti»h) 
h,e 


A  typical  term  of  B(e  ^,6^)  is 


A  D1  (e.i|/h)Dk(ef tf<h)dx 


where  A  is  piecewise  constant. 


We  have  now 


D1  (e^h^Ne^)  =  D1ejDi  (i|/^ef  )  +  z* 


where 


z*  =  .lhDk^efD1e.  +  +  e  .e^D^  - 


Az*dx  |  £  c|j|e|||2[hs  0+h2s  2°] 


Hence 
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and  so  wo  havo 


(8.41)  B(e^h,evh)  =  B(e,e<!^)  +  R 


where 


(8.42)  R<e|||u!|;2 


if  h  £  hQ(e) 


Let  us  analyze  now  the  term  BtejdJ^)  '  Obviously  supp  et^CSKD^  £)  . 
Using  theorem  5.9  and  the  remark  to  it  and  lemma  4.1,  we  can  find  w€’M(P)  , 

'Vy 

% 

supp  wca(Ph  E)  such  that 


2 

lvD(eiJvw)  l  I2l  --  cl  le,;,hi  I  1 

v€R(V)  p  h  H 1  (';)  h  H1(fl) 


l 


O/ 


% 


and  a  (V  )  <  Ko(P,  ) 
h ,  e  —  h  ,  f. 

Denote  further 


*K,«  ■  lp<*(Ph)l'pl'"(llh,c»  ■ »» 


Then  repeating  arguments  of  theorem  6.1  we  get 


(8.43) 


|B(e,i(j2e)  | 


Denote  now 
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V.  =  f A  | A C '0  ,  P€R*  } 

h ,  e  1  p  h ,  e 

Using  lemma  4  *  1  we  easily  get 

_  r\j  A/ 

(8.44)  a(Vu  )  <  Ko (V  )  <  K  a (£>  ) 

n,e  —  h,e  —  o  h ,  f 


and 


(8.45)  l 

P^,( 


!2  <  C  n2(A) 


A€P, 


h  ,e 


In  fact  if  on  w  the  bilinear  form  has  constant  coefficients  and  w  0  =  0 

P  p 

then  lemma  7.5  shows  that  ||ln  III  can  be  estimated  from  above  by  the  suitable 
turn  of  the  error  indicators.  In  the  other  cases  we  first  use  lemma-7.3  and 
lemma  7.5  and  the  same  argument  as  used  in  theorem  7.8.  Therefore  we  get 


|B(e,<f;2e)  |  <  Cf  n2  (A)  ] 1/2  |  |  ifrje  |  | 


ACJ 


h,e 


Now 


\^\ 


|  <cIU2e||  lC[||^e||2  + 

H1(n)  h  h1 (S'"  ) 

n  ,  G 


h ,  e  h  ,e 


I  l^e|  Chu-s”' 5<h),1 

h,e 

because  tj>,  ~  1  on  RJ"  and  ip,  =  0  on  fi  -  S'"  we  get 

n  n ,  c  h  h ,  c 


(8.46)  |||.i>2e|||fi  <  C(|  |^e"2 


,2  m2 


he I  I  i  +  I !^hel  I  i  6(h)  > 

H  (S'"  )  h  H1  (S'"  (h)-S'"  .) 

hjC  iijC  h  j 
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Consider 

(8.47) 

Now 

(8.48) 

and 

(8.49) 

(8.50) 

For  x€S 


now 


h  ^Hl(Sh'r(h)-Sh">f)  ^  f  (h)-S^"  c)  +  ^  ^H1(S"'  ^<h)-S"u  J 


h , e  h  ,e 


h,e  h,e 


I  K  I  ^Ho(s„,  <5(h)_s„, 
n,e  h,e 


^(S'"  6(h)-s'"  ) 

h ,e  h,e 


6(h)-sl"  )  i=l 

h,e  h,c 


H  A  (SMf  ^  V‘"/-SIM  'i  -i=i  h  ^  ct», 


H  (s’"  '  V“/-S"'  ) 
h,e  h,e 


+  |  |*  D^fce)  ||  ...  ] 

^  ^  U°/C>'l  (h)  Cll»  \ 


*T(S'M  viiy^gn,  ) 
h,  e  h  ,e 


<li,hG)  ^H°(S|;’^(h)-S^  )  <  (<<,he)^H0(S"’6(h)-S"’  ) 

n ,  e  n,£  h , e  h,e 


<  I  l^hel  I  1  C/.N 

h  H^S'"  5(h)-S"'  ) 
h ,  z  h,  e 


|  I  (D1^ ,  H .  e|  I 


S(h,  <a.-9llvll1|U  S(h).s,„ 

h,c  n,r  h,e  h ,  e 


in  ^  (h )  o'n 
ti.e  h,c 


we  have 
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dist(x,S2-S’"  *  (h))  <  h° 
h ,  e 


and  so 


IKell 


h  1  *  o 


H°^Sh'  c(h)-Sh'  c)  ^(S"'  5  (h>  _S"»  ) 


h,e  h.e2 


o(h) 

Jh,f. 

So  we  have  from  (8.43),  (8.49)  and  (8.50), 


,2  ,2 


h  ^H1(S"'  <5(h)-s"»  )  ^(S"'  6(h)-S"'  ) 

h,e  h,e  n  h,e  h,e' 


Recalling  (8.46),  (8.47),  we  get 


and  therefore 


1/2 


(8.51)  i  C  o(Dh>e)  nmaJ  |*he|  1^  <  CE1/2(„ «y |  |  V 1 1 


<  Ce 


1/2 


C  l  n2<4) 

H.  . 


1/2 


Vs  N  ,  <ccl/2||e||  ||V||  , 

H  (Q)  H  (fi)  H  (fi) 


where  we  used  the  fact  that  the  mesh  is  equilibrated  in  the  sense  of  the  assumption 
c)  above  and  that  the  estimator  is  an  upper  and  lower  one. 


Returning  to  (8.41) 


ipw* 
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ivii2i 

HX(Q) 


<  B(*. 


Ae) 


Cr1/2( 1 


i6' >H1(Q)  + 


I  I 

h  (n) 


v 


»  i  * 


and  so 


l  Vll' i 

Hx(n) 


<  Ce 


I  1 

Hx(n) 


.  e. 


B(le,i|i,e)  1  Ce  |  |e|  |21 


and  so  by  (8.40) 

IIWIlL  I  Ct||e||2 

h,£  H  (Si) 

and  the  result  follows. 

The  theorem  8.8  obviously  shows  that  the  error  estimator  is  asymptotically 
correct.  This  property  is  achieved  for  any  y  >  0  .  Let  us  discuss  now  the 
selection  of  y  .  Obviously  when  the  finite  element  solution  is  identically  zero 
then  the  term  in  the  error  indicator  associated  with  the  jumps  on  the  edges  will 
disappear  and  the  error  indicator  consists  only  of  the  "volume"  integral 

,  2  o 

y(diam  A)  £  | |g.-g  | | 

e=l  1  1  H°(A) 

2 

Assume  now  that  we  have  a  uniform  mesh  in  IR  and  the  exact  solution  is  periodic 


ui  =  Isln  IT  xi  Sln  IT  x2]ci 

r  .  2n  2 11 

u0  =  sin  v-  x,  sin  —  ~  x0]C„ 

2  hi  h  2  2 

Then  it  is  easy  to  see  that  the  finite  element  solution  is  identically  zero.  This 
leads  to  the  choice 


(8.51) 


Y 


1_  l+3u 

!lT  (X+2p)2+A2 


As  said  above,  the  choice  of  y  is  rather  arbitrary  and  (8.47)  is  one  of  the 
possibilities . 

We  have  assumed  many  very  particular  properties  of  the  meshes  and  solutions, 
in  the  theorem  8.8.  The  problem  arises  whether  these  conditions  can  be  satisfied. 
Practically  we  create  the  meshes  in  an  adaptive  mode.  The  experience  has  shown 
that  the  meshes  which  are  adaptively  constructed  have  roughly  these  properties, 
and  that  the  effectivity  index,  pj— j-|  seems  to  converge  to  1  quicker  than  the 
theorem  8.8  supports. 

2 

For  uniform  meshes  in  IR  and  smooth  solutions  it  can  be  shown  that 


HI e  I  lift  =  fc'(l+0(E2)) 


(see  [11])  and  adaptively  created  meshes  seem  to  lead  to  the  same  behaviour. 
In  the  next  section  we  will  discuss  some  concrete  illustrative  examples 


pertinent  to  these  questions. 
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9.  THE  FEARS  PROGRAM 


Based  on  the  theory  explained  above,  the  program  FEARS  (Finite  Element  Adap¬ 
tive  Research  Solver)  was  developed.  FEARS  is  a  fully  adaptive  program  solving 
a  system  of  two  elliptic  equations  and  produces  the  error  estimation  (in 
various  norms)  together  with  the  numerical  solution  of  the  given  partial  dif¬ 
ferential  equation.  The  admissible  domain  is  a  union  of  curvilinear  rectangles. 

The  adaptive  approach  is  based  on  equilibration  of  the  error  indicators.  Ihe 
description  and  experimentation  with  FEARS  will  be  reported  elsewhere. 

In  this  paper  we  are  using  FEARS  as  an  illustration  of  the  developed  theory. 

We  will  discuss  here  two  examples.  In  both,  we  are  concerned  with  the  (plane 
strain)  elasticity  problem.  We  assume  that  E  =  1  (E  is  the  Young's  elasticity 
modulus)  and  v  =  .3  (v  is  the  Poission  ratio). 

Example  1.  The  elasticity  problem  on  the  square  with  displacements  pre¬ 
scribed  on  the  boundary.  The  data  are  shown  in  Figure  9.1.  It  is  easy  to  see  that 
the  solution  belongs  to  the  space  r (fi)  (c  >  0  ,  arbitrary). 


i  ji 
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Example  2.  The  elasticity  problem  on  the  rectangle  with  mixed  boundary  con¬ 
ditions.  The  data  are  shown  in  Figure  9.2.  The  solution  is  of  the  "stamp"  type 


u  =  0 

v  =  -i 

u  =0 

V  =  1 

1  1 - — 1  1 

u  =  0 
v  =  -l 

u  =  0 

V  =  1 

1 - ! - I 

_ - _ 1 

1 - - ' - f 

1 _ l 

1 — , — * 

1 

pill  i 

i 

« 

i  1  i 

- 1 - 

Figure  9.2.  The  data  of  example  2. 

and  the  singularities  are  of  the  type  described  in  [12]  .  Solution  belong  to 
3  /  2—  e 

the  space  H  (Jl)  ,  e  >  0  arbitrary. 

In  both  cases  the  exact  solution  is  not  known,  nevertheless  by  now  elaborate 
computations  we  estimated  with  sufficient  accuracy  the  exact  energy  of  the  solu¬ 
tion.  This  gives  the  possibility  to  compute  the  (exact)  energy  norm  of  the  error 
and  compare  it  with  the  estimator. 

Example  1.  Because  of  the  obvious  symmetries  of  the  solution,  we  can  compute 
the  solution  only  on  the  quarter  of  the  original  square  applying  boundary  condi¬ 


tions  shown  in  Figure  9.3. 


Figure  9.3.  The  boundary  conditions  for  example  1. 

FEARS  constructs  adaptively  the  meshes  by  equilibrating  the  error  indicators. 

Figure  9 . 4a ,b ,c ,d ,e , f , g  show  the  sequence  of  constructed  meshes.  We  see  that 
the  sequence  of  meshes  satisfies  the  assumptions  made  in  section  8. 


Figure  9.4.  The  sequence  of  adaptively  constructed  meshes  for  example  1. 


r 
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Figure  9.5  shows  the  dependence  of  the  energy  norm  of  the  error  on  the  number 

of  cne  elements  X  .  The  norm  at  the  error  is  measured  in  percent  of  the  energy 

2—  £ 

norm  of  the  solution  (|||u|||)  .  Because  the  solution  belongs  to  H  (ft)  , 
the  rate  of  convergence  is  N  2  (or  more  precisely  N  2  )  for  the  uniform  mesh 


as  for  the  adaptive  one.  The  rate  N  2  is  the  maximal  possible  rate  because  of 
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We  aisosee  that  the  quality  of  the  error  estimator  is  better  for  the  adaptively 
constructed  meshes  than  for  the  uniform  ones.  This  is  likely  the  consequence 
of  the  equilibration  of  the  error  indicators  which  is  essential  in  our  theory. 

Example  2.  Because  of  the  obvious  symmetries  of  the  solution,  we  can 
restrict  the  computation  to  the  domain  on  a  boundary  condition  shown  in  Figure  9.7. 


Figure  9.8a,  b,  c,  d,  e,  f,  g,  h,  i  shows  the  sequence  of  the  meshes  constructed 
by  FEARS.  Once  more  we  see  that  the  assumption  about  the  mesh  made  in  section  8 
is  essentially  satisfied. 


Figure  9.8.  The  sequence  of  adaptively  constructed  meshes  for  example  2. 

Figure  9.9  shows  the  behavior  of  the  energy  norm  of  the  error.  Because  the 

n  /  <i  ^  i  <1 

solution  belongs  to  H  '  (r  >  0  arbitrary)  and  u  H  the  rate  of  convergence 

_i. 

of  the  uniform  mesh  is  N  u  .  This  is  in  complete  agreement  with  the  data  shown 

-h 

in  Figure  9.9.  The  adaptive  mesh  gives  the  rate  of  convergence  N  which  is  the 
maximal  possible  rate  for  the  smooth  solution.  We  see  that  the  adaptive  mesh 


Figure  9.9.  The  energy  norm  of  the  error  —  Example  2. 

removes  the  influence  of  the  singularities  on  the  rate  of  convergence.  We  see 
also  very  clearly  that  using  a  uniform  mesh  we  practically  can  never  achieve  an 
accuracy  of  5%. 

Figure  9.10  shows  the  behavior  of  the  effectivity  index  for  the  Example  2. 


Once  more  we  see  that  the  effectivity  index  has  practically  acceptable  value  when 
the  accuracy  of  the  solution  is  in  the  range  of  5-10%.  In  addition  the  rate  of 
convergence  of  the  effectivity  index  seems  to  be  twice  as  high  as  that  of  the 
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